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Enclosure  1 


Multiple  Antenna  Communication  with  Time  Reversal  Mirror  (MIMO-TRM) 

Pre-Processing 


Abstract 

We  have  pioneered  the  study  of  time  reversal  (TR)  pre-processing  in  radio  wireless 
communication.  Our  results  have  shown  that  TR  is  beneficial  for  use  in  a  rich  scattering 
environment  with  many  propagation  paths,  thus  positioning  the  technique  as  a  low 
complexity  solution  for.  high  data  rate  communication  via  large  delay  spread  channels. 
Methods  to  improve  the  performance  of  TR  have  been  studied,  and  comparisons  with 
other  precoding  schemes  have  been  made.  In  particular,  a  TR  system  prototype  has  been 
designed,  and  is  being  built  for  carrying  out  fijture  experiments. 

Introduction 

Communication  via  a  large  delay  spread  channel  at  high  data  rate  has  proven  to  be  a 
challenging  task,  primarily  due  to  the  substantial  amount  of  intersymbol  interference  (ISI) 
incurred  by  the  long  channel  impulse  response  (CIR).  Time  reversal  (TR)  is  a  simple  pre¬ 
processing  technique  that  makes  such  communication  feasible  without  the  need  for  any 
complex  receivers.  In  a  TR  communication  system,  the  receiver  first  sends  a  training 
sequence  to  the  transmitter  for  channel  estimation.  The  transmitter  convolves  the  message 
with  the  time-reversed  CIR,  and  sends  it  to  the  receiver.  Three  main  benefits  result  from 
this  simple  transmission  scheme. 

•  Spatial  focusing  -  The  spatial  profile  of  the  received  signal’s  power  peaks  at  the 
receiver  and  decays  rapidly  away  from  the  receiver.  This  means  that  co-channel 
interference  in  a  multi-cell  system  is  low,  and  the  probability  of  interception  by  an 
unintended  receiver  is  small.  As  a  result,  efficient  frequency  reuse  can  be  realized, 
and  secure  communication  is  attained. 

•  Temporal  focusing  -  Analogous  to  traditional  matched  filtering  at  the  receiver, 
the  power  of  all  multipaths  in  the  CIR  can  be  collected  at  one  particular  time 
instant.  Thus,  the  receiver  can  detect  the  transmitted  message  at  this  time  instant 
with  maximum  signal-to-noise  ratio  (SNR). 

•  Channel  hardening  -  The  time  variation  of  the  effective  CIR  seen  by  the  receiver 
is  significantly  less  than  that  of  the  original  channel.  In  other  words,  a  high 
frequency  diversity  gain  is  achieved,  and  the  adverse  effect  of  fading  is  reduced. 

Although  TR  has  been  studied  extensively,  and  applied  successfully,  in  the  ultrasound 
and  underwater  acoustic  arenas,  especially  in  a  rich  scattering  environment  with  ample 
multipaths,  its  suitability  for  wireless  radio  has  yet  to  be  determined.  The  core  of  our 
research  activities  was  to  investigate  the  feasibility  of  applying  TR  to  high  data  rate 
wireless  radio  communication.  The  research  results  are  summarized  in  the  following 
section.  For  detailed  explanations,  please  refer  to  the  attached  publications. 


Summary  of  Research  Results 

1.  M.  Emami,  J.  Hansen,  A.  Kim,  G.  Papanicolaou,  A.  Paulraj,  D.  Cheung,  and  C. 
Prettie,  “Predicted  Time  Reversal  Performance  in  Wireless  Communications  Using 
Channel  Measurements,”  to  appear  in  IEEE  Communications  Letters. 

This  work  represents  the  first  effort  in  estimating  the  performance  of  TR  in  wireless  radio 
using  real  channel  measurements.  Data  were  taken  in  an  indoor  environment  in  the  2  to  8 
GHz  frequency  band  by  Intel  Corporation.  Our  calculations  show  that  the  3  main  benefits 
of  TR,  namely  spatial  focusing,  temporal  focusing  and  channel  hardening,  can  all  be 
obtained.  These  results  support  the  theory  that  TR  is  an  equally  attractive  technology  for 
wireless  radio  communication  as  for  ultrasound  and  underwater  acoustic  communication. 

2.  J.  Hansen,  D.  Baum,  A.  Paulraj,  “Design  Approach  for  a  Time  Reversal  Test  Bed  for 
Radio  Channels,”  Special  Session  on  M3MO  Prototyping,  12th  European  Signal 
Processing  Conference,  Sept.  2004. 

This  work  describes  the  design  of  a  prototype  system  for  carrying  out  TR  experiments  in 
wireless  radio.  The  design  is  heavily  based  on  an  off-the-shelf  multiple-input-multiple- 
output  (MIMO)  channel  sounder,  in  order  to  keep  implementation  and  development  costs 
low.  A  company  has  tentatively  agreed  to  build  this  prototype,  pending  final  approval  of 
contract  by  Stanford  University  and  the  company.  The  availability  of  this  prototype  will 
allow  us  to  demonstrate  that  the  predicted  benefits  of  TR  can  actually  be  attained  in 
practical  environments. 

3.  T.  Strohmer,  M.  Emami,  J.  Hansen,  G.  Papanicolaou,  A.  Paulraj,  “Application  of 
Time  Reversal  with  MMSE  Equalizer  to  UWB  Communications,”  Global 
Telecommunications  Conference,  Dec.  2004. 

This  work  examines  the  performance  gain  of  adding  a  minimum  mean  square  error 
(MMSE)  linear  equalizer  at  the  receiver  of  a  TR  system.  Performance  improvement  is 
obtained  because  ISI  still  exists  in  the  received  signal,  even  though  the  signal  is  already 
time  compressed.  Moreover,  for  fixed  equalizer  complexity,  the  concept  of  rate-backoff 
is  introduced.  The  idea  is  to  reduce  the  transmission  rate,  and  hence  the  ISI,  to  a  point 
such  that  a  certain  performance  level  is  maintained. 

4.  M.  Emami,  M.  Vu,  J.  Hansen,  A.  Paulraj,  G.  Papanicolaou,  “Matched  Filtering  with 
Rate  Back-off  for  Low  Complexity  Communications  in  Very  Large  Delay  Spread 
Channels,”  submitted  to  Asilomar  Conference  on  Signals,  Systems,  and  Computers, 
Nov.  2004.  (Extended  Abstract) 

This  work  derives  the  optimum  MMSE  linear  prefilter  for  an  ISI  channel,  subject  to  the 
constraint  of  a  simple,  one-tap  equalizer  at  the  receiver.  The  performance  of  this 
transmission  strategy  is  compared  to  that  of  TR. 


5.  C.  Oestges,  J.  Hansen,  M.  Emami,  A.  Paulraj,  G.  Papanicolaou,  “Time  Reversal 
Techniques  for  Broadband  Wireless  Communications,”  European  Microwave  Week, 
Oct.  2004.  (Invited  Paper) 

This  paper  gives  an  overview  of  the  TR  technology  as  applied  to  wireless  radio 
communication.  It  contains  both  a  tutorial  introduction  and  recent  research  results, 
including  results  discussed  in  the  above  contributions. 

Conclusions  and  Future  Work 

We  have  initiated  the  investigation  on  applying  TR  to  wireless  radio  communication.  Our 
results  have  shown  that  TR  is  a  promising  technology  to  pursue,  with  the  capability  of 
achieving  high  data  rate  at  low  complexity  under  hostile  environment  with  large  amount 
of  ISI.  Future  research  will  continue  to  examine  TR  in  details.  Questions  such  as  how 
inaccurate  channel  information  affects  performance,  and  how  spatial  focusing  can  be 
quantitatively  measured,  will  be  addressed.  Comparisons  with  nonlinear  precoding 
techniques  will  also  be  made.  Moreover,  we  expect  to  perform  TR  experiments  in 
practical  environments  once  our  prototype  system  is  available.  On  a  broader  perspective, 
we  seek  to  better  understand  the  fundamental  tradeoff  between  spectral  efficiency,  power 
consumption,  performance  and  complexity  when  communicating  via  large  delay  spread 
channels  at  high  data  rate. 
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ABSTRACT 

In  time  reversal,  a  message  signal  to  be  transmitted  is  first 
convolved  with  the  time-reversed  channel  impulse  response 
and  then  sent  to  the  receiver.  The  “time-reversed”  waves 
which  propagate  into  the  channel  retrace  their  former  paths; 
they  eventually  lead  to  a  focus  of  power  in  space  and  time  at 
the  receiver.  Time  reversal  systems  can  use  antenna  arrays 
in  order  to  enhance  the  foqusing,  but  they  do  not  necessarily 
need  to.  The  channel’s  degrees  of  freedom  are  rather  excited 
by  using  broadband  signals,  and  the  focusing  can  be  achieved 
with  single-antenna  links. 

To  the  authors’  best  knowledge,  there  has  been  no 
demonstration  of  time  reversal  in  wireless  radio  yet.  This  pa¬ 
per  describes  the  broad  design  principles  of  a  system  that  is 
designed  to  demonstrate  time  reversal  in  wireless  radio  chan¬ 
nels.  The  system  block  diagram,  some  critical  specifications 
and  sources  of  error  are  described.  Experimental  data  will  be 
presented,  if  available  in  time,  at  the  workshop. 

1.  INTRODUCTION 

In  a  practical  system  employing  time  reversal  (TR),  an  in- 
tented  receiver  sends  a  training  sequence  to  an  intended 
transmitter  which  is  equipped  with  one  or  potentially  more 
than  one  transmit  antenna.  The  transmitter  time-reverses  the 
estimated  channel  impulse  response  (CIR)  and  convolves  this 
signal  with  any  message  signal  that  is  then  sent  to  the  re¬ 
ceiver. 

In  quasi-static,  reciprocal  channels  this  simple  precod¬ 
ing  scheme  yields  a  concentration  of  power  at  only  the  in¬ 
tended  receiver  at  a  particular  time.  The  spatial  and  temporal 
focusing  that  can  be  achieved  by  TR  has  been  successfully 
demonstrated  in  ultra-sound  by  Fink  [1 , 2]  and  in  underwater 
acoustics  [2, 3,  4,  5]. 

In  wireless  communications,  no  demonstration  is  known 
to  the  authors.  In  this  paper,  we  propose  to  convert  a  com- 
merically  available  broadband  MIMO  channel  sounder  into 
a  device  that  can  demonstrate  the  leverage  of  TR  in  low 
Doppler  (quasi  time-invariant)  channels  with  about  one  mi¬ 
crosecond  delay  spread.  Our  strategy  is  to  keep  the  upgrade 
as  simple  as  possible.  We  rely  on  the  MIMO  sounder’s  capa¬ 
bility  for  channel  estimation,  and  add  the  features  to  transfer 
the  CIR  estimates  back  to  the  transmitter  via  an  external  link 
and  to  then  transmit  a  precoded  message.  Hence,  we  replace 
the  pair  of  transceivers  usually  considered  in  TR  experiments 
by  a  single  transmitter  and  a  single  receiver  which  are  con¬ 
nected  by  some  feedback  loop.  The  TR  sounder  is  designed 

Jan  Hansen’s  work  was  supported  by  the  grant  HA3499/1-1  of  the 
Deutsche  Forschungsgemeinschaft  (DFG). 


such  that  this  feedback  loop  can  be  chosen  from  a  variety  of 
publicly  available,  wired  or  wireless,  links. 

The  rationals  for  this  strategy  are  twofold.  For  one,  the 
modification  of  this  type  of  existing  and  verified  system  min¬ 
imizes  the  amount  of  additional  hardware  effort  needed  and 
the  risk  of  system  performance  failure.  Secondly,  the  usage 
of  the  same  link  for  channel  estimation  and  TR  eliminates 
any  problems  with  the  reciprocity  assumptions  that  possibly 
occur  due  to  differences  in  transmit  and  receiver  RF  chains 
of  a  single  transceiver. 

The  remainder  of  the  paper  is  organized  as  follows.  In  the 
next  section,  we  give  a  brief  introduction  to  TR.  In  the  third 
section,  we  describe  the  system’s  design  principles.  We  pro¬ 
ceed  to  investigate  critial  specifications  and  sources  of  error 
in  the  fourth,  and  conclude  in  the  fifth  section. 

2.  TIME  REVERSAL 

The  transmitter  uses  the  time-reversed  complex  conjugate 
of  the  CIR  as  the  transmit  prefilter.  We  denote  the  CIR  by 
h( ro,  t),  where  ro  is  the  receiver  location  and  T  is  the  delay 
variable.  Applying  /t*{r0,-T)  as  the  prefilter,  the  effective 
channel  to  any  location  r  is  thus  given  by  the  time-reversed 
field 

s(r,T)  =  h*(r0)— T)®h(r,T)  (1) 

where  “<8>”  denotes  convolution. 

A  convolution  with  a  time-reversed  signal  is  equivalent 
to  a  correlation.  We  see  from  (1)  that  the  focusing  relies  on 
the  decorrelation  of  CIRs  in  the  delay  and  the  spatial  domain. 
A  CIR  can  be  considered  a  random  code  sequence  that  is  as¬ 
signed  to  a  transmit-receive  pair.  The  auto-  and  crosscorrela¬ 
tion  properties  of  this  sequence  are  given  by  nature.  In  con¬ 
trast  to  CDMA,  they  come  without  any  additional  bandwidth 
spreading.  The  decorrelation  of  die  CIRs  is  obviously  best  in 
a  rich  scattering  environment.  In  addition,  high  bandwidths 
unleash  the  channel’s  degrees  of  freedom  and  aid  decorrela¬ 
tion.  Physically,  one  can  consider  the  ubiquitous  locations  of 
the  scatterers  to  form  a  large  virtual  aperture  which  enables 
focusing  even  with  a  single  antenna  down  to  the  diffraction 
limit.  A  simple  but  effective  measure  for  the  focusing  capa¬ 
bility  of  a  TR  channel  is  the  delay-spread  bandwidth  product 
which  roughly  gives  the  number  of  taps  that  the  CIRs  has. 
For  an  initial  evaluation  of  the  benefits  of  TR  we  refer  to  [6]. 

3.  SOUNDER  ARCHITECTURE 

The  commercially  available  MIMO  sounding  unit  consists  of 
the  following  parts.  The  transmitter  contains  a  signal  gener¬ 
ating  unit  (SGTx)  and  a  transmitter  RF  module  (RFTx),  the 


Figure  1 :  A  schematic  of  the  TR  sounder. 


receiver  contains  a  receiver  RF  module  (RFRx)  and  a  data 
storage  unit  (DSRx).  The  conversion  of  this  sounder  into  a 
TR  sounder  is  done  by  adding  components  which  we  call 
TR  signal  processing  unit  (TRSPU)  and  feedback  loop  (FL). 
Fig.  1  shows  the  architecture  of  the  TR  sounder. 

The  proposed  TR  scheme  operates  in  three  steps;  the 
main  function  of  the  TRSPU  is  the  flexible  definition  and 
coordination  of  these  steps.  In  a  first  step,  called  channel  es¬ 
timation  step  (CES),  the  transmitter  transmits  a  training  se¬ 
quence  that  the  receiver  uses  to  estimate  the  channel.  This 
estimate  is  time-reversed  in  the  second  step,  convolved  with 
the  transmit  sequences,  and  sent  back  to  the  transmitter  via 
the  FL.  Then,  in  a  third  step  called  time  reversal  step  (TRS), 
the  transmitter  transmits  the  time-reversed  sequence  into  the 
channel;  the  receiver  now  sees  a  signal  that  is  focused  in 
space  and  time.  The  time  focusing  can  be  demonstrated  by 
the  compressed  CIR  at  the  intended  receive  antenna;  the  spa¬ 
tial  focusing  appears  because  apart  from  this  one  antenna,  all 
other  ones  will  detect  signals  that  look  like  noise. 

In  the  next  subsections  we  describe  the  specifications  of 
the  TR  sounder  in  more  detail. 


3.1  MIMO  Sounding  Unit  and  Antennas 

3.  J.  1  Bandwidth ,  Carrier  Frequency,  and  Output  Power 

The  maximum  bandwidth  of  the  currently  available  system 
is  120  MHz.  For  a  delay  spread  in  the  range  of  up  to  1  psy 
we  can  hence  expect  a  time  bandwidth  product  of  up  to  100. 
Under  normal  conditions,  wc  expect  delay-spreads  which  are 
a  factor  of  10  lower  than  that;  hence,  die  systems  works  with 
a  time  bandwidth  product  of  about  10  which  still  suffices  to 
achieve  a  reasonable  benefit  with  TR 

The  carrier  frequency  of  the  system  is  in  the  ISM  band 
at  2.4  GHz.  Unfortunately,  the  available  transmsission  band¬ 
width  is  limited  to  83  MHz.  According  to  FCC  regulations 
on  RF  devices  [7],  the  maximum  allowed  EIRP  is  4  W,  the 
maximum  transmit  power  is  1  W.  The  additional  6  dB  can  be 
achieved  by  a  suitable  choice  of  transmit  antenna. 

A  second  suitable  choice  would  be  the  ISM  band  at  5.8 
GHz  with  an  available  bandwidth  of  150  MHz.  However,  as 
the  free  space  attenuation  is  higher  and  thus  reduces  trans¬ 
mission  range,  we  decided  to  perform  the  demonstration  at 
the  lower  carrier  frequency. 

The  output  power  of  the  system  is  limited  to  1  W. 


3.1.2  Signaling 

The  sounding  signals,  and  in  particular  their  power  and 
power  spectral  density  distribution,  are  designed  in  fre¬ 
quency  domain.  The  signals  used  for  channel  estimation 
have  constant  maximum  amplitude  over  frequency  and  their 
phase  is  designed  to  minimize  the  crest  factor  and  hence 
maximize  average  output  power.  Since  all  signal  processing 
is  performed  in  frequency  domain,  instead  of  estimating  and 
time-reversing  the  CIR  the  transfer  function  (TF)  is  complex 
conjugated  in  frequency  domain. 

3.1.3  Antennas  and  Switching 

In  our  system,  the  antenna  elements  in  an  array  of  the  sounder 
are  switched.  In  this  case  the  number  of  channels  is  essen¬ 
tially  limited  by  the  number  of  available  switches  and  an¬ 
tennas  and  the  maximum  measurement  time  for  which  the 
quasi-static  channel  assumption  still  holds.  For  the  TR  ex¬ 
periments,  arrays  with  4  and  8  antennas  will  be  employed. 

Since  TR  requires  a  large  angular  spread,  omnidirectional 
antennas  are  most  advantageous.  However,  the  range  that  is 
reqiured  to  be  large  in  order  to  obtain  a  large  delay  spread 
would  benefit  from  high  gain  antennas.  The  compromise  be¬ 
tween  these  two  constraints  are  basestation  type  transmit  an¬ 
tennas  (very  narrow  vertical  and  wide  horizontal  beamwidth, 
gain  about  12dBi)  and  slightly  directional  antennas  (very 
wide  horizontal  and  wide  vertical  beamwidth,  gain  of  6  to 
10  dB). 

3.1.4  Synchronization 

The  synchronization  is  done  by  using  very  stable  reference 
clocks  at  10  MHz  (e.g.  Rubidium  standards).  The  clock  sig¬ 
nal  guarantees  longterm  stability.  In  each  RF  unit,  the  carrier 
signal  is  generated  independently  by  low-noise  PLLs. 

3.2  Feedback  Loop 

In  between  the  CES  and  the  TRS,  the  estimated  channel  data 
has  to  be  transferred  from  the  receiver  to  the  transmitter.  The 
transferral  time  must  be  much  shorter  than  the  coherence 
time  of  the  channel.  The  data  rate  of  the  link  becomes  hence 
a  critical  parameter.  A  single  complex  CIR  with  about  1000 
taps  and  8  bit  resolution  has  a  size  of  about  16  kbit.  Assum¬ 
ing  a  time  invariant  channel  with  Doppler  frequency  of  1  Hz, 
a  realistic  transmission  time  should  not  be  longer  than  0.2  s; 
consequently,  the  FL  needs  a  data  rate  of  at  least  80  kbit/s  per 
CIR  that  is  transferred. 

The  following  options  have  been  evaluated. 

3. 2 . 1  Microwave  Link 

A  microwave  link  is  advantageous  in  terms  of  data  rate, 
which  is  in  the  order  of  Mbit/s,  and  in  terms  of  flexibility.  Mi¬ 
crowave  links  can  be  established  whereeverthe  experiment  is 
to  be  performed.  However,  TR  experiments  are  likely  to  be 
conducted  in  environments  with  large  delay  spreads.  Here, 
the  TR  sounder  will  operate  at  maximum  EIRP  allowed  by 
the  FCC  regulations.  This  limit  holds  for  the  microwave  link, 
as  well.  In  areas  where  the  SNR  of  the  TR  experiment  is  low, 
that  of  the  microwave  FL  is  low,  too.  The  link  quality  can 
be  improved  by  using  directional  antennas  mounted  on  high 
poles,  but  this  requires  difficult  adjustments  and  large  efforts 
in  the  preparation  of  the  experiment 


3.2.2  Cellular  Systems 

Currently  available  cellular  systems  such  as  GPRS  offer  great 
flexibility  without  the  need  for  heavy  experimental  equip¬ 
ment.  Their  main  drawback  is  their  comparably  low  data  rate 
(order  of  100  kbit/s  for  GPRS)  and,  in  particular,  the  unreli¬ 
able  transmission  delay.  Since  the  data  rate  offered  does  not 
provide  much  overhead  compared  to  the  data  rate  required, 
the  application  of  current  cellular  services  for  the  FL  appears 
risky,  in  particular  if  the  Doppler  spread  of  the  channel  may 
at  times  exceed  1  Hz. 

3.2.3  Wired  Network 

The  wired  network  offers  tremendous  data  rates;  the  main 
drawback  is  the  inflexibility  of  its  access.  However,  most  TR 
experiments  will  be  conducted  in  urban  areas,  where  wired 
access  to  the  internet  can  be  expected  to  be  in  range.  A  wire¬ 
less  network  system  (WLAN)  such  as  one  based  on  the  IEEE 

802. 1  lb  standard  still  offers  data  rates  in  the  order  of  Mbits 
and  can  be  employed  to  bridge  the  gap  between  the  experi¬ 
mental  device  to  the  next  wired  network  connection. 

The  final  decision  about  the  best  link  depends  on  the 
availability  of  particular  services,  on  local  network  traffic, 
and  the  Doppler  spread  of  the  channel.  The  sounder  itself  is 
equipped  with  RJ45  connectors  and  can  be  hooked  up  to  any 
of  the  above-mentioned  systems. 

3.3  The  TRSPU 

The  role  of  the  TRSPU  is  to  coordinate  the  three  steps  on 
which  the  operation  of  the  sounder  is  based.  Since  it  also 
performs  the  nessecary  signal  processing,  it  is  located  in  be¬ 
tween  the  RF  and  the  data  storage  unit  of  the  receiver. 

3. 3. 1  The  Different  TR  Modes 

The  tasks  of  a  TR  sounder  are  at  least  threefold:  1)  It  must 
demonstrate  spatial  and  temporal  focusing  when  TR  exper¬ 
iments  are  performed  with  a)  a  single,  or  b)  multiple  trans¬ 
mit  antennas.  2)  It  must  be  flexible  to  allow  further  research 
about  benefits  of  TR,  i.e.  for  multiuser  (MU)  communica¬ 
tions.  3)  It  should  be  backwards  compatible.  Since  it  is  based 
on  commerical  MIMO  technology,  M1MO  channel  measure¬ 
ments  should  still  be  possible  with  a  TR  sounder.  It  turns  out 
that  4  different  TR  modes  evolve  from  these  requirements: 

1.  MIMO  Measurements: 

For  MIMO  measurements,  the  TRS  is  omitted.  The 
sounder  performs  a  predefined  number  of  channel  esti¬ 
mations  and  no  TR  transmissions. 

2.  SISO-TR: 

In  the  CES,  the  channel  estimation  is  performed  «i  times 
with  a  single  transmit  and  a  single  receive  antenna.  The 
received  TFs  can  be  averaged  to  obtain  a  better  estima¬ 
tion  of  the  channel.  The  TRSPU  complex  conjugates  this 
TF,  multiplies  it  with  the  Fourier  transform  of  the  used 
transmit  sequences,  and  sends  it  via  the  FL  to  the  trans¬ 
mitter.  The  transmitter  loads  this  sequence  and  transmits 
its  Fourier-inverse  into  the  channel.  Now,  in  the  TRS,  the 
TRSPU  provides  a  second  switching  table  and  the  RFRx 
switches  through  the  antenna  array.  At  the  one  receive 
antenna  that  is  used  in  the  CES,  the  detected  CIR  is  com¬ 
pressed.  At  all  other  antennas,  a  noise-like  signal  ap¬ 
pears.  After  «2  transmissions,  the  CES  is  repeated.  This 
mode  is  the  simplest  and  best  suited  to  demonstrate  TR. 


SGTx 


TRSPU 


Figure  2:  The  lower  layer  of  the  TRSPU,  synchronized  with 
the  SGTx. 


3.  MU-SIMO-TR: 

In  this  mode,  MU  communication  with  TR  can  be 
demonstrated.  The  CES  is  performed  with  a  single  trans¬ 
mit  antenna  and  an  antenna  array  at  the  receiver.  Concep¬ 
tually,  it  works  similar  to  the  SISO-TR  mode,  except  that 
the  TFs  of  all  receive  antennas  are  estimated.  The  TF  es¬ 
timated  at  each  antenna  is  multiplied  with  a  message  that 
is  to  be  sent  to  this  antenna;  the  sum  of  all  these  products 
is  fed  into  the  FL  and  retransmitted  in  the  TR-mode.  The 
channel  acts  now  as  a  filter;  each  receive  antenna  receives 
only  its  intended  message,  and  all  the  other  messages  as 
noise. 

4.  MISO-TR: 

In  the  original  underwater  acoustic  applications,  TR  is 
performed  with  a  transmit  array.  Each  transmit  antenna 
focuses  its  signal  on  the  same  receive  antenna,  so  that 
the  focusing  is  enhanced  proportional  to  the  number  of 
transmit  antennas.  Since  in  our  application,  the  arrays 
are  switched,  the  superposition  of  the  signals  can  only 
be  done  by  software.  Also,  the  switching  requires  short¬ 
term  phase  stability  at  the  transmitter. 

This  4th  mode  is  not  only  the  most  demanding  with  re¬ 
spect  to  synchronization,  but  also  with  respect  to  the  data 
rate  on  the  FL,  since  as  many  CIRs  must  be  transferred 
as  there  are  transmit  antennas. 

3. 3.2  Double  Layered  Structure  of  the  TRSPU 

The  TRSPU  is  implemented  in  a  double-layered  structure. 
The  lower  layer  runs  synchronized  with  the  SGTx  and  is  as 
such  displayed  in  Fig.  2.  The  SGTx  receives  a  sequence  from 
the  FL,  transforms  it  into  time  domain  and  hands  it  over  to 
the  RFTx.  The  TRSPU  receives  the  demodulated  signal  from 
the  RFRx  and  computes  the  conjugate  of  the  received  CIR 
via  an  FFT  and  division  by  the  sequence  used  for  channel 
estimation.  This  layer  is  run  in  all  modes  of  the  sounder,  be 
it  MIMO  channel  estimation,  or  one  of  the  two  stages  of  the 
different  TR  modes. 

The  upper  layer  bascially  determines  the  TR  mode,  as 
described  in  Subsection  3.3.1.  Each  mode  is  defined  by  par¬ 
ticular  numbers  n\  and  «2  for  the  duration  of  the  CES  and 
the  TRS,  and  the  one  (for  the  MIMO  mode)  or  the  two  (for 
all  other  TR  modes)  corresponding  switching  tables  that  are 
provided  to  the  RFRx.  The  upper  layer  is  software-defined. 
This  allows  any  user  to  make  changes  in  the  setup;  in  partic¬ 
ular,  the  TR  modes  can  be  embedded  in  more  sophisticated 
high-level  transmission  schemes,  e.g.,  schemes  that  adapt  the 
numbers  n\  and  «2  to  particular  channel  conditions. 

An  example  of  the  SISO-TR  mode  is  shown  in  Fig.  3.  For 


Figure  3:  The  SISO  mode  as  triggered  by  the  upper  layer  of 
the  TRSPU. 


a  given  transmit  antenna  Tx(0)  and  an  intended  receive  an¬ 
tenna  Rx(0),  the  TF  H((0,  (0, 0))  is  estimated  in  the  first  step. 
In  the  second  step,  it  is  multiplied  with  die  Fourier  transform 
of  a  transmit  sequence  and  transmitted  to  the  FL.  In  the  third 
step,  the  sounder  transmits  these  combined  sequences  until  it 
again  enters  the  first  step. 

4.  CRITICAL  SPECIFICATIONS  AND  SOURCES 
OF  ERROR 

4.1  Range  and  Power 

TR  works  best  in  environments  with  large  delay  spreads,  i.e. 
over  large  distances.  Consequently,  the  demands  on  die  link 
budget  are  high.  The  critical  parameters  are  here  the  pathless 
exponent,  the  maximum  and  average  transmit  power  deter¬ 
mined  by  hardware  and  FCC,  the  bandwidth  and  the  car¬ 
rier  frequency,  the  antenna  gains,  the  correlation  gains,  the 
sounder’s  sensitivity  and  the  target  SNR.  Initial  computa¬ 
tions  result  in  a  range  of  about  600m  for  the  described  system 
with  a  target  SNR  of  30  dB  in  a  cluttered  environment  with  a 
pathloss  of  q  —  3.5;  This  range  corresponds  to  a  propagation 
time  of  2 fis  which  is  considered  sufficient  to  create  a  delay 
spread  of  up  to  1  ps. 

4*2  Feedback  Loop 

The  FL  is  die  bottleneck  of  the  sounder.  As  discussed  in  Sub¬ 
section  3.2,  the  data  rate  required  for  the  feedback  is  about 
80  kbit  times  the  number  of  CIRs  to  be  transferred.  Depend¬ 
ing  on  the  mode,  only  one  CIR,  or  as  many  as  there  are  trans¬ 
mit  antennas,  have  to  pass  the  loop. 

If  links  with  a  data  rate  in  the  order  of  Mbits  such  as 
a  wired  network  or  a  wired  network  in  conjunction  with  a 
WLAN  are  used,  the  transmission  time  itself  is  much  shorter 
than  any  residuary  Doppler  of  the  channel  within  the  obser¬ 
vation  period.  However,  other  traffic  in  the  network  can  delay 
packet  transmission.  Initial  experiments  showed  that  packets 
of  the  size  of  64  kbit  had  travel  times  of  about  130  ms  on  an 
1 1  Mbit  WLAN  and  20  ms  on  a  local  100  Mbit  LAN. 


4*3  Phase  Noise 

TR  requires  exact  phase  information  at  the  transmitter.  Inac¬ 
curacies  in  the  phases  occur  since  a)  the  channel  is  not  per- 
feedy  time  variant,  b)  the  LO  signals  derived  by  PLLs  from 
the  10  MHz  reference  have  phase  noise. 

The  time  variance  can  be  overvome  if  the  data  rate  of  the 
FL  is  large  enough  (Subsection  3.2  and  4.2), 

Phase  noise  destroys  the  coherent  superposition  of  sig¬ 
nals  transmitted  from  multiple  antennas.  Since  it  enters  both 
at  the  transmitter  and  the  receiver,  it  will  reduce  the  received 
power  by  a  factor  of  cos(A^i  -f  A where  A&,  *  =  1,2  are 
random  variables  that  describe  the  phase  fluctuations  of  the 
transmitter’s  and  the  receiver’s  oscillators.  For  phase  errors 
less  than  15°,  this  yields  a  performance  loss  of  about  10  %  in 
SNR,  which  is  tolerable.  For  switched  antenna  sounders  this 
is  less  of  a  problem  with  respect  to  TR. 

5*  CONCLUSIONS 

This  paper  introduces  time  reversal  (TR)  as  a  promising 
transmission  technique  for  wireless  communications.  We 
propose  the  design  of  a  prototype  device  that  will  allow  the 
demonstration  of  the  benefits  of  TR.  This  prototype  is  heavily 
based  on  commercially  available  systems  and  components, 
implementation  and  development  costs  are  kept  at  a  mini¬ 
mum.  The  core  of  the  device  is  a  commercially  available 
MIMO  broadband  channel  sounder.  It  can  be  connected  to 
various  types  of  publicly  available  networks  which  serve  as 
a  feedback  loop.  Its  key  component  is  the  TR  signal  pro¬ 
cessing  unit,  which  can  operate  in  several  modes  and  offers 
flexible  handling  of  the  entire  device. 
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Abstract 

We  propose  to  apply  a  technique  called  time-reversal  to  UWB  communications.  In  time-reversal  a  signal  is  precoded  such 
that  it  focuses  both  in  time  and  in  space  at  a  particular  receiver.  Spatial  focusing  reduces  interference  to  other  co-existing  systems. 
Due  to  temporal  focusing,  the  received  power  is  concentrated  within  a  few  taps  and  the  task  of  equalizer  design  becomes  much 
simpler  than  without  focusing.  Furthermore,  temporal  focusing  allows  a  large  increase  in  transmission  rate  compared  to  schemes 
that  let  the  impulse  response  ring  out  before  the  next  symbol  is  sent.  Our  paper  introduces  time-reversal,  investigates  the  benefit 
of  temporal  focusing,  and  examines  the  performance  of  an  MMSE-TR  equalizer  in  an  UWB  channel. 

Index  Terms 

Time  Reversal,  MMSE  Equalization,  Ultra- Wideband,  Intersymbol  Interference 
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I.  Introduction 

Ultra-wideband  (UWB)  has  become  a  suitable  candidate 
for  high-data  rate,  short  range  communications.  Due  to  the 
large  operation  bandwidth,  the  resolution  in  delay  domain  is 
extraordinary,  so  that  even  in  a  dense-scattering  environment, 
the  probability  of  fading  is  low.  Since  for  the  generation 
of  these  ultra-short  pulses,  no  filters  are  required,  it  was 
considered  that  UWB  devices  are  cheap  to  produce. 

Recently,  however,  several  drawbacks  have  been  noted. 
Even  though  UWB  transmitters  irradiate  very  low  power  per 
bandwidth,  they  can  potentially  interfere  with  many  systems 
that  have  frequencies  assigned  in  the  low  GHz  range  [1]. 
Furthermore,  multipath  fading  may  be  low,  but  in  order  to 
catch  about  half  of  the  energy  that  is  distributed  in  the 
entire  impulse  response,  RAKE  receivers  with  at  least  20,  but 
potentially  many  more  taps  must  be  constructed  [2],  [3];  for 
handsets,  such  a  design  is  not  low-cost. 

We  claim  that  both  these  drawbacks  can  be  overcome 
if  UWB  is  combined  with  a  transmission  scheme  which  is 
called  time  reversal  (TR).  This  scheme  has  its  origin  in  wide¬ 
band  transmission  in  under-water  acoustics  [4],  [5],  [6],  [7] 
and  ultra-sound  [8],  [7]  and  has  recently  attracted  attention 
of  wireless  communications  engineers  [9].  In  TR,  the  time- 
reversed  channel  impulse  response  (CIR)  of  any  transmit- 
receive  link  is  taken  as  a  prefilter  at  the  transmitter.  If  such 
a  time-reversed  sequence  is  irradiated  into  the  channel,  its 
components  retrace  their  former  paths  and  lead  to  a  focus 
of  power  at  the  intended  receiver  at  some  particular  time 
instant.  For  UWB  communications,  this  strategy  has  several 
advantages.  Focusing  power  means  that  the  receiver  needs  only 
very  few  taps  to  capture  a  significant  amount  of  the  power 
in  the  channel.  Essentially,  the  complex  taks  of  estimating  a 
large  number  of  channel  taps  is  removed  from  the  receiver, 
which  may  be  a  low-cost  handset,  to  the  transmitter.  In 
addition,  the  rate  of  the  system  can  be  increased.  Whereas  with 
proposed  pulse-position  modulation  schemes  the  time  between 
the  repetition  of  two  pulses  is  chosen  to  be  sufficiently  larger 
than  the  length  of  the  CIR  [10],  focusing  of  power  means 
that  the  repetition  rate  can  be  increased.  At  last,  the  spatial 
focusing  that  comes  along  with  TR  yields  lower  interference 
with  other  communication  systems. 

In  this  paper,  we  apply  TR  to  a  single  user  downlink 
scenario  of  an  UWB  channel.  In  order  to  assess  the  usefulness 
of  TR,  we  define  a  time-compression  factor  which  measures 
the  fraction  of  the  energy  in  the  non-focused  components 
of  the  channel  versus  that  in  the  main,  focused  peak.  For  a 
channel  with  iid  complex  zero  mean  circularly  symmetric  taps, 
50%  percent  of  the  entire  CIR’s  energy  can  be  captured  with  a 
single  tap  receiver.  In  real  channels,  this  ratio  is  slightly  lower, 
which  we  show  with  the  aid  of  experimental  data  from  UWB 
measurements. 

Since  the  scattered  components  are  usually  widely  dis¬ 
tributed  over  the  delay  axis,  a  receiver  that  wants  to  capture 
more  than  50%  of  die  signal  energy  in  the  effective  CIR  must 
use  some  equalization  technique.  We  propose  an  TR-MMSE 
equalizer  for  this  purpose.  We  investigate  the  performance  of 
this  equalizer  on  real  channel  data  and  show  that  it  can  operate 


at  sampling  rates  which  are  much  higher  than  the  inverse  of 
the  length  of  the  CIR  of  the  channel;  no  ringing  out  of  the 
impulse  response  is  required. 

This  paper  is  organized  as  follows:  in  the  second  section,  we 
introduce  TR  and  define  and  investigate  the  time  compression 
factor.  Also,  the  signal  model  is  described  and  the  principles 
of  the  TR-MMSE  equalizer  are  outlined.  In  the  third  section, 
we  first  investigate  the  time  compression  factor  with  the  aid 
of  real  data.  Then,  we  examine  the  performance  of  the  TR- 
MMSE  equalizer  in  term  of  bit-error-rates  (BER).  In  particular, 
we  vary  the  number  of  taps  and  the  sampling  rate.  We  conclude 
in  Section  IV. 

II.  Theory 

A.  Time  Reversal 

The  transmitter  uses  the  time  reversed  complex  conjugate 
of  the  CIR  as  the  transmit  prefilter.  We  denote  this  CIR  by 
h(vo ,  t),  where  ro  is  the  receiver  location  and  r  is  the  delay 
variable.  Applying  the  complex  conjugate  /i*(ro,—  r)  as  the 
prefilter,  the  effective  channel  to  any  location  r  is  thus  given 
by  the  time  reversed  field 

30,  t)  =  —~h*(rQ,  -t)  ®  h(r,  r)  (1) 

where  <g>  denotes  convolution  and  the  factor  normalizes 
the  transmit  power  with  the  square  root  of  Sie  channel’s 
energy,  EH  =  /  |fi(r0,r)|2dr. 

A  convolution  with  a  time-reversed  signal  is  equivalent  to 
a  correlation.  We  see  from  (1)  that  the  focusing  relies  on  the 
decorrelation  of  CIRs  in  the  delay  and  the  spatial  domain.  In 
a  rich  scattering  environment,  a  CIR  can  be  considered  as  a 
random  code  sequence  that  is  assigned  to  any  transmit-receive 
pair.  The  auto-  and  crosscorrelation  properties  of  this  sequence 
are  given  by  nature;  unlike  in  CDMA,  they  come  without 
any  additional  bandwidth  spreading.  In  UWB,  these  sequences 
are  particularly  long.  The  benefits  of  time  compression  are 
crucial,  and  the  shortening  of  the  CIR  reduces  the  complexity 
of  equalization  at  the  receiver. 

The  decorrelation  of  the  CIRs  is  obviously  best  in  a 
rich  scattering  environment.  Physically,  one  can  consider  the 
ubiquitous  locations  of  the  scatterers  to  form  a  large  virtual 
aperture  which  enables  focusing  even  with  a  single  antenna 
down  to  the  diffraction  limit.  A  simple  but  effective  measure 
for  the  focusing  capability  of  a  TR  channel  is  the  delay  spread- 
bandwidth  product  which  roughly  gives  the  number  of  taps 
that  the  CIRs  has.  For  an  initial  evaluation  of  the  benefits  of 
TR  we  refer  to  [9]. 

B.  Temporal  Focusing 

In  a  multipath  environment,  the  channel’s  energy  Eh  is 
spread  over  delay.  It  can  be  decomposed  into  a  term  Eq 
which  describes  the  fraction  of  the  energy  in  the  direct  path, 
and  a  term  E$,  describing  that  of  the  multipaths.  In  a  strong 
multipath  environment  such  as  given  for  UWB,  we  have  the 
ratio  Es/Eo  »  1.  If  the  channel  is  slowly  time- variant,  the 
quantities  Es,  Eo,  and  Eh  are  random  variables.  Their  mean 
is  denoted  by  (.). 
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In  TR,  the  CIR  is  compressed  and  a  temporal  focus  of 
energy  is  visible  in  the  center  of  the  compressed  CIR.  In  order 
to  characterize  the  amount  of  the  temporal  focusing,  we  define 
the  temporal  compression  ratio  as 

7 tr  =  (EP)/(E™)  =  jlp  -  1  (2) 

where  EqR  is  the  energy  in  the  main  peak  of  the  received 
impulse  response,  EgR  the  one  in  the  tails,  and  Err  is  the 
sum  of  the  two. 

We  can  compute  7 tr  in  terms  of  either  the  CIR  h(r)  or 
the  transfer  function  H(co).  From  (1)  we  see  that  the  transfer 
function  Htr(uj)  of  the  time-reversed  channel  is  the  squared 
absolute  value  of  that  of  The  discrete  CIR  hTR[l],  with 
a  total  of  L  non-zero  taps,  is 

hTR[l]  =  ^  J  \H(uj)\2  exp(jujl)<iw.  (3) 

From  this  we  compute  the  energy  of  the  Oth  tap  of  the  time- 
reversed  channel, 

(e™)  =  (\n™m2)  =  flip 

=  (( Es  +  E0)2)  =  (E2h)  (4) 


and  the  energy  of  the  entire  channel, 

{Ethr)  =  I h[l]  ®  h*[-l}\2)  =  I*  |tf(w)|4dw)  (5) 


where  the  frequency  domain  representation  on  the  right  hand 
side  of  the  equation  follows  from  Parseval’s  theorem.  Whether 
a  time  domain  or  a  frequency  domain  representation  is  more 
desirable  depends  on  in  which  domain  information  about  the 
channel  is  available. 

If  we  evaluate  expression  (5)  further  in  time  domain  for  L 
identical  and  independently  distributed  taps  with  equal  power, 
we  can,  derive  the  time  compression  factor  as 

04  +  2  (L 


7th  = 


fa  +  {L-  1)/?! 


i. 


(6) 


Here,  02  and  04  denote  the  second  and  the  fourth  moment  of 
the  amplitude  in  each  tap,  respectively.  For  L  —  1,  one  has 
7TjR  =  0;  with  increasing  L  it  converges  to  1 .  Hence,  in  the 
time-reversal  channel,  the  ratio  of  the  energy  of  the  zeroth  tap 
to  the  sum  of  the  energy  of  all  the  other  taps  converges  to  a 
fixed  value. 

If  /i[Z]  is  a  complex  Gaussian  random  variable  with  mean 
0  and  variance  (E0  +  Es)  =  { Eh ),  we  have  02  =  (Eh)  and 
04  =  2 (Eh)2-  Hence, 


iTR 


2  L(Eh)2  __  2 

2(Eh)2  +  {L-1)(Eh)*  1  +  i 


(7) 


For  large  L,  the  power  in  the  tails  of  the  impulse  response 
becomes  about  equal  the  power  in  the  peak.  Hence,  even 
though  ISI  can  be  significantly  suppressed,  any  system  without 
equalization  operates  at  best  at  0  dB  signal  to  ISI  plus  noise 
ratio.  For  small  L ,  the  compression  works  in  the  mean  better 
than  for  a  large  one,  even  though  the  difference  quickly  levels 


Fig.  1.  Schematic  of  TR  impulse  hm  and  truncated  TR  CIR  h^R.  In  this 
toy  example  L  —  29,  Ls  =9  and  the  number  of  effective  taps  (marked  as 
fi  lied  circles)  is  5. 


out.  But  the  probability  that  a  channel  realization  occurs  that 
does  not  focus  at  all  (i.e.,  all  channel  taps  are  real  and  have 
the  same  magnitude)  is  much  larger  than  for  long  sequences. 

C.  Signal  Model  and  TR-MMSE  Equalizer 

The  transmitter  modulates  symbols  Sk  using  a  pulse  shaping 
filter  0(r)  of  bandwidth  B  —  l/Ts .  It  then  performs  TR  on 
the  pulse  stream.  Thus: 

x(t)  =  Y,  sk<t>(t  -  kMTs)  <g>  h*{-t)  (8) 

k 

where  h(r)  is  the  (infinite  bandwidth)  CIR.  After  matched 
filtering  the  received  signal  can  be  written  as: 

y(t)  =  '£2sk<f>(t  -  kMTs)  ®h*(-t)  <g)  h(t)  <8<p*(-t)  -F  n(t) 

k 

(9) 

The  receiver  then  samples  y(t)  at  IMTS : 

y[l]  =  y(lMTs)  =  s^hTR[l  -  *]  +  n[l],  (10) 

k 

where  hrR(t)  =  h(t)  <E)  h*(— t),  h(t)  =  h(t)  0  <f>{t),  and 
fiTR[k ]  =  hxR(kMTs).  The  sampled  received  signal  is  thus 
the  convolution  of  the  symbols  with  a  downsampled-by-M 
version  of  Iitr{.)- 

Due  to  the  temporal  focusing  described  above  the  equalizer 
at  the  receiver  does  not  have  to  use  the  quite  long  (time- 
reversed)  CIR.  Instead  the  equalizer  uses  a  “shortened”  TR 
CIR  by  keeping  only  a  few  taps  which  capture  most  of  the 
energy  of  Iitr  and  still  achieve  very  good  performance.  More 
precisely,  we  construct  a  shortened  TR  CIR  h^R  of  length  L3 
as  illustrated  in  Fig.  1  by  defining 

KM  S  ,S'  OK 

for  |Z|  <  Lsj 2  and  Ls  —  max{/  :  \Iitr[1]\  >  6}.  Here  5  is 
an  a  priori  chosen  tolerance  level  that  depends  on  the  target 
SNR.  For  convenience  we  will  rename  the  indexing  of  h^R 
such  that  hj,R  —  . . . ,  h^R[Ls  —  1]]*. 
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The  symbols  are  transmitted  one  block  at  a  time  with  a 
guard  period  of  length  Ls.  The  output  of  the  channel  is: 

y[m]  =  (x  0  hTR)[m]  4-  nR[m] 

Here  x[m\  is  the  input  signal  m  =  0, . . . ,  N  —  1  and  y[m]  is 
the  output  signal  m  =  0, . . . ,  N  +  L  -  2.  In  matrix  notation 
this  can  be  written  as  Toeplitz-type  system 


y  =  H  thx  +  nR 


where  y  —  [y[ 0] . . .  y[N  4-  L  -  2]]T  and  x  =  [x[0] . . .  z[i\T  - 
1]]T.  Using  h^R  instead  of  hrR  for  the  MMSE  equalizer  we 
replace  HtR  by  the  matrix  Hf  R  of  size  (N  +  Ls  —  1)  x  N 
given  by 

f  0  0  1 


H 


s  _ 

Tfi  — 


~  1] 

0 


0 


0  ...  0  —  1]  J 


and  the  matrix  equation  for  the  MMSE  estimator  becomes 

*  =  (H^)*y. 

Assuming  the  transmitted  symbols  are  from  a  BPSK  alphabet, 
we  can  now  decode  x  by  looking  at  the  sign  of  each  element. 

We  note  that  is  a  sparse  Toeplitz  matrix,  thus  x  can 
be  computed  efficiently  by  combining  sparse  matrix  techniques 
with  fast  Toeplitz  solvers  [1 1]. 


III.  Simulations 

A,  Measurement  Data  and  Temporal  Focusing 

Measurements  were  conducted  with  a  network  analyzer  by 
Intel  Corp.  at  off-peak  hours  to  ensure  channel  stationarity. 
The  environment  is  an  office  space  (40m  x  60m)  with  many 
cubicles;  measurements  were  conducted  at  several  locations. 
They  span  a  bandwidth  of  2-8  GHz  with  3.75  MHz  frequency 
resolution.  Antennas  are  vertically  polarized.  The  measure¬ 
ments  are  corrected  to  compensate  for  the  system  components 
(including  cable,  gain  stages,  and  antennas).  The  height  of 
the  transmit  antenna  is  about  2.5m  and  that  of  the  receive 
antenna  is  lm  above  the  floor.  As  a  sample  result,  Fig.  2 
shows  the  CIR  (above)  and  the  time-reversed  CIR  (below). 
The  compression  is  clearly  visible.  Still,  this  figure  is  not 
suitable  to  tell  the  ratio  between  the  power  in  the  peak  of 
the  time-reversed  impulse  response  and  that  in  the  tails.  In 
order  to  demonstrate  this  ratio,  we  plotted  in  Fig.  3  the  time 
compression  ratio  y TR  for  a  CIR  length  between  1  and  500 
taps,  where  the  entire  length  of  the  CIR  is  fixed  at  l/3.75MHz. 
The  curve  was  actually  computed  in  frequency  domain;  in 
order  to  suppress  fluctuations  due  to  the  randomness  of  the 
samples,  a  sliding  window  was  implemented,  so  that  7 Tr  was 
actually  averaged  over  all  available  frequency  samples. 
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Fig.  2.  Impulse  Response  of  the  channel  (above)  and  of  the  time-reversed 
channel  (below) 


Fig.  3.  In  the  TR  channel,  the  ratio  between  the  eneregy  in  the  main  tap  and 
the  energy  in  the  tail  of  the  CIR  converges  for  independent  taps  to  a  fi  xed 
ratio. 


For  a  small  number  of  taps,  the  curve  coincides  very  well 
with  the  theoretical  one.  At  about  200  taps,  the  maximum 
theoretical  value  of  1  is  exceeded;  beyond  that,  correlation 
between  the  taps  starts  to  play  a  role  and  the  power  in  the 
tails  of  the  CIR  will  be  stronger  than  that  in  the  main  peak. 


B.  Performance  of  the  Equalizer 

We  used  the  single  channel  realization  whose  magnitude 
response  at  full  bandwidth,  B,  is  shown  in  Fig.  2.  Three 
cases  were  considered.  Case  I  has  M  =  8  (see  Eqn.  10)  and 
bandwidth  B.  Case  II  has  M  —  4  and  bandwidth  B/2.  Finally 
case  HI  has  M  =  2  and  bandwidth  B/A.  The  constellation 
was  BPSK  and  no  channel  coding  was  used.  Thus  all  the 
scenarios  have  the  same  data  rate  and  all  channels  have  250 
taps.  Fig.  4  shows  the  bit  error  rate  (BER)  curves  for  these 
scenarios  when  the  receiver  only  estimates  the  20  strongest 
taps.  When  M  =  2  the  curve  floors  very  rapidly  since  not  a 
very  significant  portion  of  the  energy  of  the  CIR  is  captured  by 
the  largest  20  taps.  When  M  =  4,  the  ISI  energy  is  due  only  to 
every  fourth  tap  of  the  TR  channel.  Thus  more  of  the  energy 
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Fig.  4.  BER  vs  SNR  when  the  receiver  only  estimates  the  20  largest  taps 
of  the  channel.  Bandwidth  expansion  factor  of  8  (solid  line),  4  (dotted  line), 
and  2  (dashed  line) 


Fig.  5.  BER  vs  SNR  when  the  receiver  only  estimates  the  10  largest  taps 
of  the  channel.  Bandwidth  expansion  factor  of  8  (solid  line),  4  (dotted  line), 
and  2  (dashed  line) 


Fig.  6.  BER  vs  SNR  when  the  receiver  only  estimates  20  largest  taps  of  the 
channel  -  another  channel  realization.  Bandwidth  expansion  factor  of  8  (solid 
line),  4  (dotted  line)  and  2  (dashed  line). 


Fig.  7.  BER  vs  number  of  strongest  taps  estimated  at  the  receiver.  SNR  = 
15  dB 


of  the  effective  CIR  is  captured  by  the  strongest  20  taps.  This 
effect  improves  significantly  in  the  case  where  M  =  8.  Note 
that  the  case  M  —  8  uses  more  bandwidth.  However,  the  BER 
floors  at  a  reasonable  SNR  of  15  dB  for  UWB.  Fig.  5  shows 
the  BER  curves  for  the  same  three  cases  but  when  the  receiver 
estimates  only  10  strongest  taps  of  the  CIR.  In  this  scenario 
the  error  flooring  occurs  earlier  for  all  the  cases.  However,  the 
complexity  of  the  receiver  is  significantly  reduced. 

Fig.  6  plots  the  BER  for  a  different  CIR.  Here  the  receiver 
estimates  the  largest  20  taps  of  the  CIR.  Comparison  between 
this  figure  and  Fig.  4  shows  that  the  BER  performance  is  very 
similar  for  these  two  channel  realizations.  This  is  due  to  the 
large  diversity  available  in  this  ultra- wideband  channel. 

Fig.  7  shows  the  BER  vs  number  of  strongest  CIR  taps  that 
the  receiver  estimates  for  equalization  at  15  dB  SNR.  This  plot 


shows  that  the  BER  reduces  approximately  exponentially  with 
increasing  number  of  taps.  Fig.  8  shows  the  BER  vs  number 
of  strongest  CIR  taps  that  the  receiver  estimates  at  15  dB 
SNR  when  the  transmitter  does  not  use  TR.  Figs.  7  and  8 
demonstrate  the  temporal  focusing  ability  of  TR.  The  receiver 
requires  to  estimate  much  less  taps  with  TR  than  without  it  in 
order  to  achieve  the  same  BER  performance. 

Note  that  the  plots  are  for  the  strongest  taps.  These  taps 
may  not  be  next  to  one  another.  There  may  be  zeros  between 
any  two  consecutive  strongest  taps.  Thus  the  total  length  in 
samples  of  the  estimated  channel  is  in  general  larger  than  the 
number  of  strongest  taps  estimated.  However,  the  complexity 
of  equalization  is  approximately  the  same  as  if  the  length  of 
the  estimated  channel  is  equal  to  the  strongest  number  of  taps 
estimated. 
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Fig.  8.  BER  vs  number  of  strongest  taps  estimated  at  the  receiver  with 
regular  transmission.  SNR  -  15  dB 


IV.  Conclusion 

Time  reversal  significantly  reduces  the  number  of  taps  that 
the  receiver  needs  to  estimate  in  order  to  achieve  a  certain 
target  BER.  If  the  CIR  has  iid  complex  Gaussian  taps  the 
ratio  of  signal  to  ISI  power  approaches  unity  as  the  length  of 
the  CIR  increases  without  bound.  We  show  that  this  ratio  is 
very  close  to  unity  when  the  CIR  is  more  than  about  10  taps 
long. 

In  order  to  mitigate  the  ISI  one  can  fix  the  bandwidth  and 
increase  symbol  spacing  by  an  integer  factor  of  M.  The  larger 
M,  the  less  taps  of  the  CIR  the  receiver  needs  to  estimate  in 
order  to  attain  a  reasonable  BER.  However,  larger  M  means 
more  wasted  bandwidth.  With  no  channel  coding  and  by  only 
estimating  the  20  strongest  taps  out  of  approximately  250 
channel  taps,  the  receiver  can  reach  10" 3  BER  at  around  12  dB 
SNR  with  MMSE  equalization  and  M  —  8.  In  order  to  achieve 
such  a  performance  without  TR,  the  receiver  needs  to  estimate 
almost  all  the  taps  of  the  CIR. 
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Extended  Abstract 

We  study  the  possibility  to  transmit  data  over  channels  with 
large  delay  spreads  under  the  constraint  of  a  very  simple  re¬ 
ceiver  which  has  only  one  tap.  Such  a  scheme  is  of  interest 
when  a  cost-efficient  way  to  transmit  potentially  large  data  rates 
is  sought.  We  investigate  the  performance  of  the  optimal  pre¬ 
filter  for  this  scheme,  and  compare  it  to  a  simplified,  so-called 
time-reversal,  prefilter  which  has  very  low  complexity.  Time 
reversal  filters  have  in  connection  with  wireless  transmission 
over  large  bandwidths  their  origin  in  underwater  acoustics  and 
ultrasound  [1].  Their  application  to  wireless  communications 
was  recently  demonstrated  on  UWB  data  [2].  At  low  SNR,  die 
optimal  prefilter  and  the  TR  prefilter  are  equivalent.  At  high 
SNR,  the  time  reversal  prefilter  achieves  a  performance  that  is 
independent  from  the  delay  spread  of  die  channel  and  hence  the 
same  for  any  bandwidth. 

In  this  contribution  we  derive  analytically  the  performance 
of  the  discussed  system.  We  discuss  the  impact  of  multiple 
transmit  antennas  and  of  a  transmission  scheme  that  uses  rate 
back-off,  i.e.,  it  transmits  at  a  rate  lower  than  the  inverse  of  the 
bandwidth.  If  bandwidth  is  abundant,  such  as  it  is  the  case  in 
Ultra-wideband  communications  [3],  our  proposed  scheme  can 
achieve  a  reliable  quality  of  service  at  high  data  rates  with  a 
single  tap  receiver.  We  mention  that  our  setup  corresponds  to 
the  downlink  of  an  access  point  to  a  handheld  In  principle, 
our  analysis  holds  equivalently  for  a  simple  transmitter  with  a 
receive  equalizer  with  a  large  number  of  taps. 


A.  Prefiltering  for  Long  Delay  Spread  Channels 

A  block  diagram  of  the  system  under  study  is  given  in  Fig.  1 . 
The  transmit  symbols  x  are  fed  into  equalizers  and  then  trans¬ 
mitted  through  the  channels  described  by  their  channel  impulse 
responses  (CIRs)  hi9i  =  1. . . . ,  Mr ,  and  Mp  is  the  number  of 
transmit  antennas.  The  noise  is  denoted  by  n,  and  the  receiver 
has  a  single  tap  whose  gain  is  denoted  by  (3.  Since  with  very 
large  delay  spread  the  transmit  prefilter  is  not  likely  to  eliminate 
inter  signal  interference  (ISI)  completely,  we  allow  that  the  sys¬ 
tem  does  not  transmit  at  full  rate  (inverse  bandwidth),  but  at  a 
integer  fraction  of  this  rate,  which  we  call  D.  The  input-output 


Fig.  1 .  A  block  diagram  of  a  MISO  system  with  a  preequalizer  at  the  trans¬ 
mitter  and  a  single  tap  receiver. 


relation  of  this  system  can  be  written  as 

Mt 

y[D][k ]  =  - n+Mk]  (1) 

*=1  l 

where  [A;]  is  the  output  symbol  at  time  k.  The  rate  back-off 
factor  D  =  1 , . . . .  L  is  equivalent  to  upsampling  the  signal  be¬ 
fore  transmission  by  a  factor  of  D,  and  L  is  the  channel  length. 

We  derive  the  optimal  prefilter  solving  the  minimization 
problem 

g p,N  =  argmaxg  ^_E[[|  xt/JHg  +  (3n\k\  -  x\k  -  A]  ||2] 
s.t.  E  ||  g*g  ||=  1.  (2) 

The  product  Hg  is  the  convolution  of  the  prefilter  and  the  chan¬ 
nel,  the  latter  stacked  into  a  toeplitz  matrix,  and  N  is  the  length 
of  the  filter.  The  vector  g  is  given  by 

g=[g[0]T  ...  g[JV-l]T]T  (3) 

and 

g[«]  =  [»[«]  ■■■  9mMT  1  =  0,. ..,N- 1.  (4) 

A  is  the  delay  of  the  equalizer  and  is  approximately 
The  difference  between  this  formulation  and  the  one  where  the 
equalizer  is  placed  at  the  receiver  is  the  transmit  power  con¬ 
straint.  The  above  problem  is  convex  in  g  and  ,3  and  can  be 
solved  analytically  using  the  method  of  Lagrange  multipliers: 

fr.jv  =  I  (hhH  +  HHeA  (5) 

f32=  H*eA 


(6) 
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We  furthermore  define  the  matched-filter  bound  SNR  Pmfb  as 


PMFB  =  - ~2~ - •  (9) 

Combining  (8)  and  (9)  we  have  that  the  effective  SNR  of  a  TR 
communication  system  with  Mt  transmit  antennas  and  a  rate- 
back  off  of  D  is 


Peff,D  =  PMFB 


I+PMFB^^ 


(10) 


In  our  paper  we  will  present  a  formula  for  the  ratio  of  the  use¬ 
ful  power  against  the  ISI  power,  E[Ptr]/ E[Q^]  for  a  chan¬ 
nel  with  exponentially  decaying  power  delay  profile  and  delay 
spread  crT.  Using  this  result,  we  derive  the  following  limiting 
cases: 


Fig.  2.  The  effective  channel  as  generated  by  TR  (dashed  line)  and  the  MMSE 
fi  her  (coni  line) 


Eq.  (5)  defines  a  filter  with  an  MMSE-type  solution  with  a 
finite  number  of  taps.  We  notice  that  the  optimized  value  for  (3 
bascially  ensures  that  the  power  constraint  of  die  filter  is  valid. 
In  the  transmit-receiver  chain  of  Fig.  1,  f3  has  no  effect  on  the 
SNR  at  the  receiver;  it  ensures,  however,  that  the  transmitter 
can  take  into  account  the  SNR  to  optimize  the  preequalizer. 

We  introduce  a  simplified  filter  which  we  call  TR-prefilter 
and  which  basically  is  a  matched  prefilter;  it  follows  from  the 
optimal  filter  (5)  for  low  SNR, 


g TR 


lim  g  =  h* 

SNR-> 0 

s.t.  II  h  ||=1, 


(7) 


or 

0T7l,»[fc]  =  y/lK  [-fc], 

with  7  =  ^£^  £f=-11  |/ii[i]|2)  ,  i.e.,  the  filter  equals  the 

time-reversed  CIR  of  the  channel,  and  the  receiver  receives  the 
autocorrelation  function  of  the  channel.  The  power  in  its  central 
peak  is  maximized,  but  the  ISI  suppression  is  worse. 

The  CIR  of  a  pre-equalized  channel  with  originally  50  taps 
and  constant  power  delay  profile  is  shown  in  Fig.  2;  for  the 
computation  of  die  MMSE  filter,  p  ~  Ptx/&2  =  5  dB  was 
assumed.  The  dashed  line  is  obtained  from  the  MMSE  filter,  the 
continuous  one  from  the  TR  prefilter.  We  see  that  the  prefilters 
compress  the  CIRs  significantly  around  a  strong  peak.  The  ISI 
of  the  impulse  response  obtained  with  TR  is  higher  than  that  of 
the  MMSE  filter,  but  its  peak  is  slightly  higher,  too.  For  a  large 
number  of  taps,  the  peak  is  essentially  non-fading  which  yields 
reliable  system  performance. 


B.  Performance  Evaluation  of  the  TR  filter 

For  performance  evaluation,  we  denote  by  Pq  and  the 
instant  useful  power  and  the  ISI  power,  respectively,  and  define 
the  effective  SNR  pe/f  as 

PTxE[Po]  _  E[P0] 

Peff’D  PrxE{QM]+a 2  J5[QP1]  +  I-  W 


Peff,D 

Peff,D 


Pmfb  for  °r  — *  0 

DMtpmfb  o 

- — —  for  a. 

Pmfb  +  DMt 


(H) 

oo  (12) 


Eq.  (11)  tells  that  a  TR  system  is  lossless  should  the  delay 
spread  of  the  channel  be  very  small.  Eq.  (12)  can  be  split  into  a 
low-SNR  and  a  high- SNR  regime  as 


peff,D  =  Pmfb  for  pmfb  low.  aT  large  (13) 
Pe.ff,D  =  DMt  for  pmfb,  crT  large.  (14) 


Eq.  (13)  shows  that  TR  is  optimal  at  low  SNR;  this  behavior  is 
expected  for  a  transmit  matched  filter.  Numerical  results  pre¬ 
sented  in  the  final  paper  will  show,  however,  that  the  low-SNR 
region  ends  at  about  0  dB  which  not  sufficient  for  reliable  com¬ 
munications.  Eq.  (14)  demonstrates  that  in  an  environment  with 
a  large  delay  spread  aT,  the  performance  of  a  TR  system  at 
high  input  SNR  saturates  at  a  value  independent  from  the  delay 
spread  of  the  channel.  The  performance  increases  with  3  dB  per 
added  transmit  antenna,  and  with  another  3  dB  for  each  rate- 
back  off  of  a  factor  2  from  the  maximum  available  bandwidth. 
We  will  also  include  simulations  that  show  that  the  central  tap 
is  non-fading. 

If  bandwidth  is  abundant,  we  can  hence  design  a  commu¬ 
nication  system  with  a  very  simple  receiver  that  has  only  one 
tap  and  does  not  need  to  perform  channel  estimation.  This  tap 
is  non-fading  and  the  effective  SNR  is  not  a  function  of  the 
bandwidth  of  the  system.  Varation  of  transmission  rate  or  us¬ 
ing  some  particular  number  of  transmit  antennas  guarantees  any 
reasonable  quality  of  service.  In  our  final  paper  we  will  discuss 
the  derivations  and  results  presented  here  more  thoroughly,  give 
some  numerical  examples  which  highlight  some  of  the  analyti¬ 
cal  results,  and  point  out  future  applications  of  the  investigated 
system. 
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Abstract — Channels  with  large  delay  spreads,  which  typi¬ 
cally  arise  when  large  bandwidths  are  used,  are  considered  as 
disadvantageous  for  wireless  communications.  However,  despite 
their  challenges,  these  channels  also  offer  opportunities.  The 
application  of  time  reversal  techniques  may  be  one  possibility 
to  exploit  these  advantages  and  to  make  communications  in 
channels  with  large  delay  spreads  more  feasible.  This  paper 
describes  advantages  and  draw  backs  of  time  reversal  techniques 
in  broadband  wireless  communications,  and  outlines  concepts, 
goals  and  potential  improvements.  The  properties  of  time  re¬ 
versal  in  various  simulated  and  measured  wireless  channels  are 
demonstrated,  and  system  design  is  discussed. 

I.  Introduction 

While  there  is  little  doubt  that  future  communication  sys¬ 
tems  will  have  to  accomodate  larger  data  rates,  die  natu¬ 
ral  solution  consisting  of  increasing  the  bandwidth  brings 
a  number  of  open  challenges.  In  this  matter,  time  reversal 
(TR)  techniques  could  play  a  role  in  the  design  of  broadband 
wireless  communications.  In  TR,  channel  state  information 
(CSI)  at  one  or  possibly  more  transmit  antennas  is  used  to 
precode  transmitted  symbols  with  the  time  reversed  version 
of  their  respective  channel  impulse  responses  (CIRs).  These 
“time  reversed”  waves  propagate  in  the  channel,  retrace  their 
former  paths  and  eventually  lead  to  a  focus  of  power  in  space 
and  time  at  the  receiver  [l}-~{3].  hi  rich  scattering  media,  this 
space-time  compression  can  be  very  strong,  and  at  the  focal 
point,  the  effective  channel  obtained  through  TR  is  hardened. 

Spatial  focusing  [4]  means  that  the  spatial  profile  of  the 
power  peaks  at  the  intended  receiver  and  decays  rapidly  away 
from  the  receiver.  Temporal  focusing  means  that  the  CIR  at  the 
receiver  has  a  very  short  effective  length.  Channel  hardening 
is  the  phenomenon  of  tightening  of  the  distribution  function 
of  the  effective  CIR’s  power. 

Using  TR  for  wireless  communications  in  channels  with 
large  delay  spreads  presents  therefore  three  merits. 

•  The  temporal  focusing  compresses  the  channel’s  energy 
into  one  single  peak  with  very  low  sidelobes.  The  task  of 
equalizing  the  channel  is  hence  significantly  simplified. 

♦  The  spatial  focusing  can  result  in  much  lower  intercept 
probability  and  co-channel  interference  in  a  multi-cell 
system.  By  selectively  focusing  the  energy  in  both  space 
and  time  at  a  target  point,  TR  ensures  that  intercept 


receivers  will  have  difficulty  detecting  or  decoding  the 
intended  signal. 

•  Through  the  channel  hardening,  TR  provides  high  diver¬ 
sity  gain.  The  statistics  of  the  time  reversed  channel  are 
different  from  the  actual  channel.  Specifically  the  time 
reversed  channel  has  a  much  smaller  variance  than  the 
physical  channel  itself. 

Time  reversal  is  well  known  in  acoustics  and  has  lead  to  re¬ 
markable  applications  in  underwater  sound  [5]— {13]  and  ultra¬ 
sound  [14]-{16].  The  extension  of  TR  techniques  to  wireless 
communications  has  yet  to  be  precisely  investigated,  although 
the  idea  of  exploiting  scattering  has  recently  been  proposed 
in  [17] — [20].  From  a  communication  engineer’s  viewpoint,  TR 
shifts  the  system’s  matched  filter  from  the  receiver  side  into 
the  channel.  The  complexity  of  the  transmitter  increases;  the 
receiver  captures  a  preequalized  CIR  which  is,  in  the  case  of 
large  delay  spread  channels,  much  simpler  to  equalize  than  the 
original  one.  Research  on  TR  touches  a  number  of  different 
research  areas.  Seen  as  a  preequalizer,  TR  is  known  as  a 
transmit  matched  filter  [21]  or  pre-rake  [22]  which  has  been 
studied  in  connection  with  CDMA  but  not  at  very  large  delay 
spreads,  where  spatial  and  temporal  focusing  comes  into  effect. 
TR  has  to  be  compared  to  optimal  schemes  (e.g.  [23],  [24]). 
Complexity  will  become  a  major  issue  since  delay  spreads  of 
100  or  more  taps  can  make  many  known  optimal  schemes 
infeasable.  The  spectral  efficiency  is  also  of  great  interest. 
TR  does  not  use  the  capacity  achieving  power  allocation  [25], 
but  in  fading  channels,  an  increased  bandwidth  yields  channel 
hardening  which  could  increase  outage  capacity.  Finally,  many 
practical  issues  need  to  be  addressed.  At  which  bandwidth 
does  it  become  advantageous  to  precode  transmitted  messages 
in  order  to  shorten  file  channel,  and  which  precoding  method 
should  be  applied?  At  which  size  of  a  network  does  it 
become  useful  to  shift  complexity  from  several  receivers  to 
one  centralized  transmitter?  We  believe  that  the  recent  interest 
in  transmission  over  very  large  bandwidths  such  as  provided 
by  the  standard  IEEE  802. 15  [26]  is  a  good  incentive  to  rethink 
many  of  the  long  known  methods  that  were  invented  to  deal 
with  channels  of  low  and  moderate  delay  spread.  We  propose 
TR  as  a  suitable  starting  point  for  these  investigations,  as  it 
clearly  points  out  the  potential  gains  that  high-delay  spread 
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channels  may  offer  despite  their  draw  backs. 

Here,  we  shall  present  an  initial  study  of  TR  applied  to  radio 
wireless  systems.  In  Section  II,  we  briefly  introduce  TR,  to¬ 
gether  with  the  formalism  used  in  this  paper.  In  Section  III,  we 
compare  the  rate  of  a  TR  system  to  one  that  uses  waterfilling 
and  we  simulate  the  cumulative  distribution  function  of  the  rate 
of  such  a  system  in  a  fading  channel.  Section  IV  deals  with 
the  characterization  of  time-reversed  random  fields  in  wireless 
channels.  To  this  end,  we  first  define  suitable  metrics  in  order 
to  characterize  temporal  and  spatial  focusing.  Then  we  present 
measured  results  of  TR  space-time  focusing  in  an  indoor  Ultra 
Wide  Band  (UWB)  radio  channel  in  the  2-8  GHz  band  [28]. 
To  provide  more  insight  in  the  space-time  focusing  results,  we 
finally  perform  simulations  for  broadband  systems  (bandwidth 
up  to  100  MHz)  using  two  very  different  channel  models  at 
2.5  GHz  [27]. 

In  Section  V,  we  investigate  the  applicability  of  TR  to 
wireless  data  transmissions.  We  examine  how  TR  performs 
compared  to  optimized  linear  precoding,  and  we  compute  bit 
error  rates  of  both  a  TR  and  an  optimized  system  [29].  We 
describe  in  Section  VI  how  commercially  available  MIMO 
channel  sounders  can  be  converted  into  sounders  that  can  be 
used  for  real-world  TR-channel  measurements,  and  we  finally 
conclude  in  Section  VII. 

II.  Time  Reversal  Formulation 

Consider  a  transmission  between  transmit  point  P  and 
receive  point  Q.  The  CIR  is  denoted  as  /i(r,P  — >  Q).  It 
is  important  to  note  that  the  usual  symmetry  properties  imply 
that  /i(r,  P^Q)  =  /i(r,  Q P). 

The  so-called  CIR  is  actually  the  convolution  of  the  infinite- 
bandwidth  physical  channel  response  and  the  filter  impulse 
response  of  bandwidth  B ,  e.g.  a  Nyquist  filter  with  given  roll¬ 
off  factor. 

Depending  on  the  scattering  channel  and  the  bandwidth, 
the  CIR  results  in  a  temporal  spreading  of  the  initial  pulse. 
Scatterers  indeed  create  multipath  mechanisms  which,  in  turn, 
cause  echoes  to  arrive  at  the  receiver  with  different  delays. 
However,  the  resolvability  of  the  different  delays  depends 
upon  the  ratio  of  the  inverse  of  the  bandwidth  to  the  physical 
channel  spread  (i.e.  the  interval  between  successive  delays). 
The  smaller  \/B  is  relative  to  the  channel  delay  spread, 
the  larger  is  the  number  of  resolved  paths.  As  an  important 
conclusion  of  the  above  discussion,  it  can  be  expected  that  the 
B  x  7rms  product  will  play  a  significant  role  in  TR  focusing 
results. 

Using  the  proposed  formalism,  the  received  signal  at  any 
point  R  for  a  pulse  emitted  from  R0  is  h(r,Ro  — ►  R). 
At  a  particular  point  T,  which  we  define  as  the  transponder, 
h(r,  Ro  — »  T)  is  captured.  If  the  transponder  sends  back  the 
time  reversed  version  of  the  captured  signal,  i.e.  h(— r,  R0  — > 
T),  then  at  any  point  R,  the  received  signal  can  be  expressed 
as 

s(r,  R)  =  h{r,  T  R)  *  h(-r,  R0  T)  (1) 

where  ★  denotes  the  convolution  product.  Note  that  (1)  as¬ 
sumes  a  perfect  estimation  of  h(r,  R0  — >  T).  In  practical 
settings,  noise  and  interference  considerations  will  cause  this 


estimation  to  be  biased.  It  should  also  be  remembered  that 
s(t,  R)  is  related  to  the  bandwidth  B. 

From  now  on,  we  define  the  point  R  =  R0  as  the  focal,  the 
target  point  or  simply  the  receiver.  Furthermore,  when  there  is 
no  ambiguity  about  the  transmitter  position,  we  also  drop  the 
symbol  T  and  the  arrow  in  the  notation  of  the  CIR. 

Based  on  (1)  and  on  the  symmetry  properties  mentioned 
above,  the  signal  received  at  Ro  is 

s(r,Ro)  =  A(r,Ro)*/i(-r,Ro)  (2) 

As  a  consequence  of  (2),  the  time  reversal  operation  causes 
the  received  signal  at  Ro  to  be  focused  in  both  time  and  space 
through  constructive  interference,  meaning  that  all  multipath 
signals  add  coherently  at  Ro,  and  incoherently  elsewhere. 

Models  usually  compute  the  channel  transfer  function 
H(uj,  Ro)  over  the  system  bandwidth  (where  uj  is  the  angular 
frequency).  It  is  the  Fourier  transform  of  h(r,  R0).  Relation¬ 
ships  (1)  and  (2)  are  easily  written  in  the  frequency  domain, 
since  the  time  reversal  operation  corresponds  to  a  complex 
conjugation  in  frequency  (denoted  by  the  superscript  *): 

S(u,  Ro)  -  H(lj,  Ro )H*{u,  Ro)  (3) 

According  to  (3),  we  observe  that  the  time  reversal  operation 
is  equivalent  to  a  perfect  channel  matching,  i.e.,  the  channel 
acts  as  the  system’s  matched  filter.  According  to  (2),  we  see 
that  s(r,  Rq)  is  the  autocorrelation  function  of  the  channel,  so 
that  there  is  a  second  interpretation  for  TR.  We  may  consider 
the  CIR  as  a  randomly  generated  code  sequence  which  is 
given  by  nature.  Then,  a  TR  communication  system  operates 
similar  to  a  broadcasting  CDMA  system  with  a  pre-RAKE 
[22].  But  the  code  sequence  is  not  a  spreading  sequence 
in  the  sense  that  it  spreads  transmission  bandwidth  without 
increasing  information  rate  for  each  single  user.  In  TR,  the 
fully  spread  bandwidth  can  be  used  for  transmission  to  each 
user.  However,  since  the  sequence  does  not  have  carefully 
designed  correlation  properties  but  is  at  best  random,  the 
TR  system  cannot  erase  intersignal  interference  (ISI)  and 
multiuser  interference  completely. 

III.  Capacity  and  Outage  in  TR  Schemes 

We  compare  in  this  section  the  rate  of  a  TR  system  to  one 
that  uses  waterfilling  and  simulate  the  cumulative  distribution 
functions  of  the  rates  of  these  systems  in  a  fading  channel. 
TR  enforces  a  particular  power  allocation  scheme,  P(uj)  — 
\H(u>)\2,  which  is  in  general  not  capacity  achieving.  The 
power  allocation  strategy  inherent  to  TR  is  to  enhance  the 
strong  components  of  the  channel,  and  to  weaken  those  which 
are  less  strong.  At  low  SNR,  this  strategy  is  related  to  water¬ 
filling,  so  that  we  can  expect  a  close-to-optimal  performance. 
At  high  SNR,  where  a  constant  power  allocation  scheme  is 
optimal,  TR  is  not.  We  demonstrate  this  behavior  in  Fig.  1 
in  which  the  mean  rate  of  a  TR  communication  system  that 
operates  in  a  channel  with  50  i.i.d.  Gaussian  distributed  taps 
is  depicted  (dashed  line).  The  rate  equals  at  low  SNR  the 
capacity  of  the  channel  (continuous  line),  and  deviates  more 
and  more  at  higher  SNR.  Also  for  comparison,  the  capacity 
of  the  AWGN  channel  is  shown.  It  is  largest  in  this  picture 
since  all  curves  are  normalized  to  a  bandwidth  of  1  Hz  [25]. 
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IV.  Time  Reversal  Channel  Characterization 


A.  Characterization  Metrics 

The  space-time  received  signal  $(r,R)  is  converted  into 
performance  metrics  by  considering 

•  the  energy  of  s(t,  R)  at  any  point  R  in  space,  at  a  given 
time  ro,  i.e. 

t?d(R)  =  |s(r0,R)|2  (4) 


with  ro  such  that  |a(r0,Ro)|  =  maxr{|s(r,Ro)|}. 

•  the  RMS  delay  spread  of  s(t,  R)  (on  a  realization  basis), 


At(R)  - 


f(r 


1)2|s(r,  R)|2dr 


/|s(r,R)|2dr 


(5) 


Fig.  I.  The  average  rate  of  a  TR  communication  system  reaches  at  low  SNR 
the  capacity  of  the  frequency  selective  channel. 


where  rm  is  the  average  delay  defined  as 

TVn(R)  =  J  t\s(t,  R)|2dr/  J  |s{r,R)|2dr. 


Cumulative  Distribution 


Note  that  the  above  definitions  assume  that  the  received 
space-time  signal  is  sampled  correctly  in  time  so  that  the 
maximal  amplitude  can  be  captured.  An  alternative  definition 
of  7?( R)  would  be  to  consider  the  peak  energy,  irrespective  of 
the  time  delay: 

Vm(R)  =  max{|s(r,  R)|2}.  (6) 

T 

Both  tjd( R)  or  tjm (R)  and  Ar(R)  are  random  spatial  func¬ 
tions,  which  can  be  characterized  by  their  first-order  moments. 

1)  Spatial  Focusing:  The  spatial  focusing  around  the  focal 
point  is  described  by  two  parameters. 

The  asymptotic  space-focusing  gain  (SFG)  is  given  by 

pd  =  iFf  »u>(R'o)/»to(H.).  (7) 

|R— Ro|— ►oo 


Fig.  2.  In  channels  with  a  larger  number  of  taps,  the  outage  probability 
decreases. 


In  fading  channels,  outage  capacity  can  be  more  relevant 
that  ergodic  capacity.  The  cumulative  distribution  function 
(cdf)  of  the  rate  of  a  TR  system  at  10  dB  SNR  in  channels 
with  10,  50,  and  100  i.i.d.  Rayleigh  taps  is  shown  in  Fig.  2. 
The  cdf  of  the  TR  system  is  shown  as  the  blue  curves,  those 
of  a  scheme  that  uses  water  filling  in  black.  The  rate  of  the 
TR  system  has  for  a  particular  number  of  taps  a  cdf  which 
is  very  similar  to  that  of  the  optimal  scheme  but  shifted  by  a 
fixed  amount  of  less  than  a  quarter  bps/Hz  to  the  left.  A  TR 
system  in  a  channel  with  50  taps  has  at  an  outage  probability 
of  20%  about  the  same  rate  as  an  optimally  precoded  system 
that  uses  water  filling  in  a  channel  with  10  taps.  The  increase 
from  50  to  100  taps  does  not  improve  the  rate  at  this  outage 
probability  much. 


It  is  the  ratio  between  the  energy  at  Ro  to  the  energy  at  long 
distance  from  Ro.  A  large  value  of  this  ratio  indicates  better 
space  focusing.  Note  that  pM  is  defined  similarly  with  respect 
to  7?Af  (R),  and  that  pd  could  be  defined  at  different  time 
delays,  i.e.  other  than  r0. 

The  3-dB  contour  of  the  energy  function  77#  ( R)  or  77 a/ (R) 
can  be  considered  as  the  focusing  region.  It  is  described  by 
the  distance  in  both  range  and  cross-range  for  which  ?7d(R) 
or  t}m  (R)  remains  within  3  dB  below  the  energy  at  Ro.  The 
characteristic  parameters  Ga  and  Gx  are  therefore  defined  (e.g. 
for  r)o  such  that: 

7?d(Ro  +  Ga  ua)/7?D(Ro)  =  0.5  (8a) 

t}D(Ro  4-  Gx  ux)/t/d(Ro)  -  0.5  (8b) 

where  ua  and  u*  are  unit  vectors,  respectively  in  the 
range  and  cross-range  directions  (i.e.  in  directions  respectively 
parallel  and  orthogonal  to  the  T  -  R  axis). 
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2)  Time  Focusing:  The  time  focusing  at  the  focal  point  is 
described  by  the  RMS  delay  spread  of  s(r ,  R)  at  R  —  R0, 
denoted  as  Aro  =  At(Ro).  Note  that  this  del  ay- spread  is 
expressed  in  (5),  and  accounts  for  the  pulse  width. 

Finally,  a  time  focusing  gain  (TFG)  is  also  suitably  defined 
by  the  relative  increase  of  RMS  delay  spread  Ar(R)  at  any 
point  R  compared  to  Rq,  The  parameter  is  denoted  as 


a(R)  = 


Ar(R)  -  Arp 
A  r0 


(9) 


The  asymptotic  TFG  is  given  by  limjR_Ro|_>00  <r(R).  A  larger 
TFG  indicates  better  temporal  focusing  in  the  sense  that  the 
time  compression  at  the  focal  point  with  respect  to  any  point 
away  from  the  focal  point  becomes  larger. 

The  performance  measures  descrived  above  are  random 
variables,  depending  on  the  channel  realization.  Hence,  they 
should  be  characterized  by  their  mean  and  variance  taken  over 
all  realizations. 


B.  Measured  UWB  Channels 

In  this  first  part,  we  present  CIRs  obtained  from  TR  as  com¬ 
puted  from  real  data  measured  with  a  very  large  bandwidth. 

1)  Experimental  Set-Up:  Measurements  were  conducted  by 
Intel  Corp.  at  off-peak  hours  to  ensure  channel  stationarity. 
The  environment  is  an  office  space  (40  m  x  60  m)  with 
many  cubicles.  Measurements  span  the  bandwidth  2-8  GHz 
with  3.75  MHz  frequency  resolution.  Antennas  are  vertically 
polarized.  The  virtual  grid  on  which  the  receiver  is  moved  has 
a  distance  of  Am/4  where  Am  is  the  wave  length  of  the  mid 
frequency  of  the  measurements  (5  GHz).  The  antenna  is  moved 
to  a  different  location  with  a  precise  robotic  positioner.  At  each 
antenna  position,  the  channel  is  measured  over  a  6  GHz  band 
from  2.0  GHz  to  8.0  GHz  using  a  vector  network  analyser. 
The  measurements  are  corrected  to  compensate  for  the  system 
components  (including  cable,  gain  stages,  and  antennas).  The 
height  of  the  transmit  antenna  is  about  2.5  m  and  that  of  the 
receive  antenna  is  lm  above  the  floor. 

The  CIR  is  computed  by  taking  the  inverse  DFT  of  mea¬ 
surement  data. 

2)  Experimental  Results:  We  evaluated  t?m  for  1 1  different 
scenarios.  Hereafter,  we  only  consider  the  two  extreme  cases 
which  have  the  worst  focusing  and  yet  they  are  very  similar 
to  the  outcomes  for  the  other  scenarios.  The  first  case  is 
an  essentially  line-of-sight  (LOS)  scenario  with  a  transmitter- 
receiver  separation  of  3  m.  The  second  case  has  a  separation 
of  1 1  m  and  is  a  typical  non-LOS  (NLOS)  situation. 

A  3D  plot  of  77m  (R)  of  these  two  cases  is  shown  in  Fig.  3. 
In  this  figure,  the  square  grid  spans  a  region  of  7Am  x  7Am 
where  Am  «  6  cm.  We  see  that  spatial  focusing  works  fine 
in  both  scenarios.  In  neither  of  them,  however,  the  peak  is 
isotropic;  both  peaks  have  one  direction  in  which  they  fall  off 
faster,  and  another,  in  which  the  decay  is  slower.  The  structure 
of  the  peak  carries  some  information  about  the  geometry  of 
the  environment,  i.e.,  about  directions  which  show  faster  and 
those  which  show  slower  decorrelation  in  space. 

The  signal  power  level  is  at  least  10  dB  lower  at  a  distance 
of  7Am  than  its  value  at  file  receiver.  This  demonstrates  the 
spatial  decorrelation  very  well. 


£ 


(b)  NLOS  scenario 

Fig.  3.  Spatial  Focusing  r; m(R)-  One  shot  spatial  field  realizations  for  a) 
the  line-of-sight  b)  the  non-line-of-sight  scenario. 


In  Fig.  4  we  illustrate  the  time  compression  of  the  received 
pulse.  Displayed  are  the  magitude  of  the  CIR  for  the  LOS 
scenario  (a),  the  respective  time-reversed,  i.e.,  compressed  CIR 
(b),  the  CIR  for  the  NLOS  scenario  (c),  and  the  compressed 
one  for  this  case  (d).  The  channel’s  power  that  is  not  in 
one  (central)  tap  is  significantly  shortened.  Yet,  for  the  LOS- 
case,  comparably  strong  side  lobes  remain  visible,  since  the 
deterministic  component  causes  correlation.  The  CIR  of  the 
NLOS  channel  is  much  better  compressed  in  time.  Still, 
there  is  some  intersymbol  interference  in  the  time  compressed 
field  and  this  would  require  equalizing  the  received  signal. 
However,  the  time  compressed  channel  has  a  significantly  less 
power  in  its  side  lobes  and  is  thus  less  expensive  to  equalize. 

For  these  two  sets  of  data,  we  computed  the  temporal  and 
the  spatial  focusing  gain  which  are  displayed  in  Tab.  I.  The 
TFG  is  75%  and  30%,  respectively,  for  the  LOS  and  the  NLOS 
case.  The  CIR  with  file  stronger  deterministic  component  has 
a  higher  TFG.  The  SFG  is  21  and  42  so  that  the  SFG  in  the 
LOS  scenario  it  is  lower.  A  dominating  strong  component  in 
the  channel  yields  higher  TFG,  but  lower  spatial  orthogonality. 

To  provide  a  better  physical  understanding  of  these  observed 
results,  we  will  present  in  the  following  a  detailed  character- 
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(a)  LOS 


(b)  LOS 


TABLE  I 

TFG  AND  SFG  COMPUTED  FROM  EXPERIMENTAL  DATA. 


TFG  (eq.  (9)) 

PM  (eq.  (6),  (7)) 

LOS 

0.75 

21 

NLOS 

0.31 

42 

or  urban  canyons. 

Here  multipaths  are  caused  by  two  different  mechanisms: 
reflections  from  the  waveguide  walls  and  scattering  by  the 
medium  inhomogeneities.  Fig.  5  shows  a  sketch  of  the  waveg¬ 
uide:  two  rays  are  propagating  from  the  source  at  R0  to  the 
transponder  T.  One  ray  reaches  T  after  several  reflections 
from  the  waveguide  walls.  The  other  one  reaches  T  by 
multiple  scattering  on  the  inhomogeneities.  In  fact  the  source 
at  Rq  produces  an  entire  angular  spectrum  of  rays  over  the 
bandwidth,  and  signals  reaching  T  undergo  a  mixture  of  these 
two  multipathing  effects. 


500  lOOO  1500 

taps 


(c)  NLOS 


Fig.  4.  Magnitude  of  channel  impulse  responses  (a  and  c)  and  the  time 
compressed  impulse  response  (b  and  d) 


ization  of  TR  focusing  in  two  different  wireless  channels  at 
2.5  GHz:  a  randomly-filled  waveguide  and  a  discrete  scattering 
channel. 

C.  Random  Waveguide 

1)  Channel  Model:  We  investigate  in  this  section  guided 
wave  propagation  in  two  dimensions.  The  interior  of  a  periodic 
waveguide  is  a  weakly  scattering  medium.  Signals  propagating 
through  this  waveguide  experience  significant  multipathing 
yielding  a  large  delay-spread.  This  model  can  seen  as  an  ap¬ 
proximation  of  propagation  through  hallways  inside  buildings 


\  v 


7T 

“V1  "-VT 


'>v 

v  “;V V/; 

- 


A 


• 

f 

\J 

_ 


VI 

T 


Fig.  5.  A  sketch  of  the  random  waveguide  mode!  with  two  rays  from  the 
source  at  Ro  propagating  to  the  transponder  at  T  each  undergoing  one  of 
the  two  mechanisms  for  multipathing:  refection  from  walls  and  scattering  by 
inhomogeneities. 


The  waveguide  width  is  Lx  and  the  distance  along  the 
waveguide  between  R0  and  T  is  Lz.  We  assume  that  the 
initial  transmitted  signal  at  Ro  can  be  given  a  limited  angular 
aperture  d.  For  the  sake  of  simplicity,  the  mean  refractive  index 
is  chosen  as  unity  (note  that  changing  this  value  only  affects 
the  overall  phase  of  the  computation).  The  fluctuation  is  an 
isotropic,  Gaussian  correlated  random  function  of  space,  with 
RMS  height  ms  and  correlation  length  £. 

To  compute  the  random  space-time  field,  we  place  equi- 
spaced  phase  screens  inside  the  waveguide  separated  in  dis¬ 
tance  by  L  Simulating  propagation  along  the  waveguide  in¬ 
volves  combining  a  sequence  of  steps  from  one  screen  to 
the  next.  Each  step  involves  propagation  through  an  “empty” 
waveguide  followed  by  a  random  phase  correction. 

We  now  describe  a  single  step  of  the  phase  screen  method. 
Let  zn^i,  zn  and  zn+ 1  denote  locations  of  three  consecutive 
phase  screens.  We  assume  that  the  channel  transfer  function 
JIb(w,  zn-i  — ►  zn)  from  zn_ i  to  zn  is  known.  Because  the 
waveguide  is  periodic,  we  express  it  as  a  Fourier  series: 


oo 

H(w,zn-i  —  2n)  =  am(<jj,zn)e>2nmx/L*.  (10) 

m= — oo 
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The  channel  transfer  function  for  zn  to  zn+i  is  then  given  by 


H(uj,zn^  zrl+l)  =  ejk>i(-x'>5z/2 

OO 

J2  Tm(yJ,5z)am{u,zn)e*™x/L* 


(11) 


with  Sz  =  zn+i  —  zn  and  k  =  2n/ A  the  wavenumber.  The 
forward  propagation  operator  Tm(u),  Sz)  is  defined  as 


Tm(u,Sz)  =  exp 


jkyj  1  -  (27t m/(kLx)): 


Sz 


(12) 


The  random  phase  correction  ji{x)  is  a  path-integral  of  the 
fluctuation  fj.(x ,  z)  =  1  -  n2(x ,  z)  between  the  screens  zn  and 
zn+ i  with  n(x,  z)  denoting  the  index  of  refraction  inside  the 
waveguide. 

To  limit  the  angular  aperture  for  the  initial  transmit  at  Ro, 
we  filter  the  Fourier  modes  before  updating  the  field  to  the 
next  phase  screen.  Let  zq  denote  the  phase  screen  containing 
Ro  and  z0,x)  denote  the  initial  transmit  signal  which 
we  express  as  the  Fourier  series: 


OO 

<B(uj,zo,x)  =  ^2  ipm(uj)e?2'*rnx/L*.  (13) 

m=— oo 


From  (12)  we  determine  that  propagating  modes  are  those  for 
which 

H<tt-  <14> 

All  others  are  evanescent.  For  a  propagating  mode,  the  prop¬ 
agation  angle  is  determined  from  sin  =  2irm/(hLx).  The 
largest  value  of  m  such  that  (14)  holds  gives  the  largest 
propagation  angle  $max-  To  set  the  angular  aperture  of  the 
initial  transmit  we  restrict  the  summation  in  (14)  to  modes 
m  such  that  <  i9. 

2)  Simulation  Results :  Numerical  simulations  in  the  waveg¬ 
uide  have  been  carried  out  using  the  following  parameters: 

•  the  central  frequency  is  2.5  GHz 

•  the  length  Lz  is  100  m 

•  the  refractive  index  fluctuation  has  RMS  height  s  — 
0.05  and  correlation  length  t  —  1.0  m 

•  the  waveguide  width  Lx  is  either  12.8,  25.6  or  51.2  m 

•  the  source  angle  is  either  25°,  45°  or  65° 

•  the  bandwidth  B  is  between  10  and  100  MHz. 

The  focus  point  Ro  is  centered  in  the  waveguide.  We  set 
Ro  to  be  the  origin  for  the  x  —  z  coordinate  system.  The 
transponder  is  located  100  m  down-range  from  R0.  It  is  an 
aperture  spanning  1  m  <  x  <  2  m.  Note  that  the  aperture 
center  is  deliberately  slightly  skewed  from  x  —  0  as  to  break 
any  symmetry.  Hence,  any  observed  effect  is  by  no  means 
an  artifact  caused  by  symmetry.  Note  also  that  our  numerical 
results  have  shown  that  the  space-time  focusing  is  insensitive 
to  the  location  of  the  transponder  array  as  long  as  there  is 
sufficient  multipathing  in  the  channel. 

Figs.  6  and  7  show  r/p(x,  Zo)/r)D(xo7z0)  and  A r(x,z0) 
respectively  for  a  waveguide  of  width  25.6  m.  Averages  were 
computed  for  500  realizations  of  the  random  refractive  index. 
Results  are  shown  for  bandwidths  of  10,  20,  40,  80  and  100 
MHz. 


Fig.  6.  Energy  of  sb(to,x)  in  a  waveguide  with  width  Lx  —  25.6  m  and 
Lz  =  100  m  averaged  over  500  realizations  of  the  random  refractive  index. 
The  source  is  limited  in  angle  by  A#  =  45°  and  in  bandwidth  by  100  80, 
40,  20  and  10  MHz. 


Fig.  7.  RMS  delay  spread  trms  averaged  over  500  realizations  of  the  random 
refractive  index.  All  other  parameters  are  the  same  as  in  Fig.  6. 


In  Fig.  6  we  observe  that  quality  of  spatial  focus  about  x0  = 
0  is  nearly  the  same  over  all  bandwidths.  However,  the  energy 
away  from  the  refocus  location  decreases  as  the  bandwidth 
increases.  In  particular,  we  observe  a  5  dB  difference  between 
the  10  MHz  and  100  MHz  cases. 

In  Fig.  7  we  observe  the  quality  of  temporal  focus  about 
#o  —  0  where  At  is  the  smallest.  Because  a  larger  bandwidth 
yields  a  shorter  pulse.  At  is  smaller  for  larger  bandwidths. 
Nonetheless,  the  curves  become  smoother  as  the  bandwidth 
increases.  This  phenomenon  is  due  to  the  onset  of  statistical 
stability  manifested  from  broad  bandwidth  signals  [34]. 

To  examine  spatial  focusing  in  greater  detail,  we  show  the 
space-focusing  gain  in  Fig.  8.  The  top  plot  is  for  Lx  =  25.6 
m  with  different  values  of  The  bottom  plot  is  for  $  =  45° 
with  different  values  of  Lx.  As  we  have  already  mentioned, 
adding  bandwidth  to  the  system  reduces  the  energy  away  from 
the  refocus  location.  Fig.  8  further  demonstrates  this  since  all 
curves  show  that  SFG  increases  with  bandwidth.  The  top  plot 
in  Fig.  8  shows  that  limiting  reduces  the  SFG.  The  bottom 
plot  shows  that  there  is  little  difference  between  the  Lx  =  12.8 
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Fig.  8.  Asymptotic  space-focusing  gain  (SFG)  as  a  function  of  bandwidth 
for  different  source  angles  (top)  with  waveguide  width  Lx  —  25.6  m  and 
different  waveguide  widths  (bottom)  with  source  angle  A#  =  45°. 


bandwidth  (MHz) 


Fig.  9.  3-dB  contour  in  cross-range  Gx  as  a  function  of  bandwidth  for 
different  source  angles  (top)  with  waveguide  width  Lx  =  25.6  m  and  different 
waveguide  widths  (bottom)  with  source  angle  At?  =  45°. 


m  and  25,6  m  cases.  However  die  SFG  for  the  Lx  =  51.2 
m  case  is  much  smaller  than  the  other  two.  For  a  fixed 
widening  the  waveguide  reduces  the  diversity  of  modes  that 
propagate  in  the  channel  This  reduction,  in  turn,  reduces  the 
quality  of  focus  in  time  reversal. 

In  addition,  we  show  the  3-dB  contour  in  cross-range  Gx 
in  Fig.  9.  Similar  to  Fig.  8,  the  top  plot  is  for  Lx  —  25.6 
m  with  different  values  of  &  and  the  bottom  plot  is  for 
=  45°  with  different  values  of  Lx.  Here  we  observe 
that  cross-range  dependence  of  the  3-dB  contour  changes 
very  little  with  respect  to  bandwidth.  Because  limiting  the 
angular  aperture  of  the  source  is  a  spatial  low-pass  filter 
yielding  broader  initial  transmit  sources,  Gx  increases  as  # 
decreases.  For  different  waveguide  widths,  the  3-dB  contour 
does  not  change  significantly.  It  decreases  slightly  as  the  width 
increases  because  of  added  angular  diversity.  From  (14)  we 
understand  that  the  number  of  propagating  modes  within  a 
fixed  angular  aperture  is  set  by  A  and  Lx.  With  A  fixed  and  Lx 
increasing,  more  propagating  modes  exists  inside  the  angular 
aperture.  For  example,  with  d  ~  45°,  there  are  approximately 
75,  150  and  301  propagating  modes  for  Lx  =  12.8  m,  25.6  m, 
and  51.2  m,  respectively.  Time  reversal  methods  exploit  this 
added  spatial  diversity  to  tighten  the  refocusing. 


To  examine  temporal  focusing  more  closely,  we  show  the 
asymptotic  limit  of  the  time- focusing  gain  in  Fig.  10.  The 
top  plot  is  for  Lx  —  25.6  m  with  different  values  of  t? 
and  the  bottom  plot  is  for  ^  —  45°  with  different  values  of 
Lx.  All  curves  decrease  monotonically  with  bandwidth.  The 
main  reason  for  this  decrease  is  that  the  pulse  is  narrower 
as  bandwidth  increases.  The  overall  delay  spread  is  reduced 
and  the  refocus  region  is  not  as  pronounced.  The  TFG  in¬ 
creases  with  d  since  die  additional  diversity  allows  for  greater 
temporal  focusing  in  the  sense  that  ratio  between  A r  at  and 
away  from  xo  becomes  greater.  Similar  to  the  result  for  the 
SFG,  we  observe  that  the  TFG  for  Lx  =  12.8  m  and  25.6  m 
are  close  to  each  other  while  the  widest  waveguide  yields  a 
smaller  TFG.  As  Lx  increases  effects  manifested  by  reflections 
from  waveguide  walls  decrease  and  the  waveguide  domain 
approaches  an  unbounded  one  exhibiting  a  smaller  TFG. 

D.  Geometry-Based  Stochastic  Channel 

1)  Channel  Model:  For  discrete  scattering  media,  such 
as  typical  outdoor  wireless  environments,  the  channel  can 
be  represented  by  a  number  of  effective  scatterers  randomly 
distributed  in  space.  A  ray-based  approach  can  then  be  used  to 
describe  file  channel  as  a  sum  of  so-called  scattered  or  reflected 
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Fig.  10.  Asymptotic  time-focusing  gain  (TFG)  as  a  function  of  bandwidth 
for  different  source  angles  (top)  with  waveguide  width  Lx  —  25.6  m  and 
different  waveguide  widths  (bottom)  with  source  angle  A#  =  45°. 


contributions.  In  the  following,  the  statitical  distribution  of 
scatterers  is  based  on  the  well-known  one-ring  model  [30]. 

We  first  consider  that  only  single  scattering  occurs.  We 
assume  that  the  scatterers  are  uniformly  distributed  inside  an 
annular  region  surrounding  the  target  point,  as  illustrated  in 
Fig.  1 1 .  This  region  is  specified  by  a  minimal  radius  pm  and 
a  maximal  radius  pm-  The  angle-spread  of  the  channel  (as 
seen  from  file  target  point)  is  fixed  by  limiting  the  aperture  of 
the  annular  region  to  a  given  portion,  i.e.  specifying  minimum 
and  maximum  angles,  and  i 9m  (see  Fig.  1 1).  In  this  paper, 
we  refer  to  Atf  =  —  dm  as  the  scattering  angle-spread  or 

simply  the  angle-spread. 

We  further  simplify  the  channel  description  by  choosing 
Pm  ~  2pm  and  pM  «  D,  with  D  denoting  the  distance 
between  the  transmit  and  target  points.  Hence,  the  RMS  delay 
spread  of  the  physical  channel  (over  an  infinite  bandwidth) 
is  well  approximated  for  an  omnidirectional  distribution  of 
scatterers  (A#  =  27 r)  by 


trms 


2e 


(15) 


with  c  —  3  x  108m/s  is  the  speed  of  light.  For  a  limited  angle 
spread,  we  assume  that  —  7r-  Atf/2  and  =  7T+ Atf/2. 


Fig.  11.  Geometrical  representation  of  the  propagation  model. 


In  that  case,  the  delay-spread  can  be  written  as 
\/5pm  /  At?  +  sin  Atf 

TRMS*—  V - 2tt -  (16) 

The  channel  transfer  function  is  easily  estimated  by  means 
of  a  ray-tracing  approach  that  yields 


H{to,  Ro) 

N 

l  — l 

rfRoS^SjT  ] 


L^H<3Si  +  dstT 


7/2 


GfaSt^T)  (17) 


N 


-Er« 


exp  [-ja;(<feoS,  +  rfs,T)] 


[^RoS,  +  <fe,T 


17/2 


(18) 


The  quantities  in  (17)  and  (18)  are 

•  (j  =  27r/A,  the  angular  frequency  (A  denoting  the 
wavelength); 

•  (^(u;,  P  — ►  Q)  =  exp(— ju;dpQ)/(dpQ)7/2,  file  spread¬ 
ing  function  for  a  transmission  from  P  to  Q,  with 
7  denoting  the  path-loss  exponent  (7  =  2)  and  dpg 
denoting  the  distance  between  P  and  Q; 

•  F/,  the  scattering  coefficient  of  the  2th  scatterer;  in  this 
paper,  it  is  considered  to  be  a  complex  Gaussian  variable, 
with  given  standard  deviation; 

•  usua^  Spreading  factor  for  reflec- 
non  [32]  ,  ensuring  that  the  transfer  function  is  inversely 
proportional  to  the  total  distance,  as  outlined  by  (18); 
multiplying  this  reflection  spreading  factor  by  T*  yields 
the  so-called  scattering  cross-section. 

Note  that  (17)  and  (18)  implicitly  assume  that  the  CIR 
is  filtered  by  a  Nyquist  filter  with  roll-off  factor  equal  to 
zero,  hence  the  filter  is  a  rectangular  window  in  the  frequency 
domain.  Simulations  were  done  with  the  following  parameters: 

•  the  central  frequency  is  2.5  GHz; 

•  the  transponder  is  separated  from  the  target  point  by  a 
distance  of  104  wavelengths,  i.e.  1.2  km; 

•  when  the  delay-spread  is  not  used  as  a  variable  or 
explicitely  assigned  a  specific  value,  it  is  fixed  at  0.8/zs 
(i.e.  pm  =  1800A); 

•  different  angle  spreads,  ranging  from  7r/3  to  27t. 

•  the  scattering  coefficient  is  assumed  to  be  complex  Gaus¬ 
sian  distributed  with  standard  deviation  of  0.25. 
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2)  Omnidirectional  Scattering  Results:  Fig.  12  shows  snap¬ 
shots  of  t)d( R)  (the  target  point  is  in  the  center  of  the 
figure),  for  two  bandwidtbs:  50  and  100  MHz,  two  channel 
delay  spreads  and  omnidirectional  angular  spreading  of  the 
scatterers. 


two-(X4jim-a=10»»te  xm*“16**  B  =  lOUH/ 


Fig.  12.  Typical  one-shot  realizations  of  time  reversed  random  fi  elds  (the 
energy  function  is  expressed  m  [dB]). 

The  impact  of  both  the  bandwidth  and  the  channel  delay- 
spread  is  clearly  visible  (the  scales  are  kept  constant  from 
graph  to  graph).  For  die  smallest  B  x  trms  product,  some  areas 
receive  the  same  level  of  energy  as  the  target  point.  For  the 
largest  B  x  trms  product,  the  space  focusing  gain  is  about  15 
to  20  dB.  Also,  it  seems  that  the  instantaneous  energy  function 
is  oscillating  as  a  function  of  the  distance  to  the  focal  point 

Fig.  12  shows  that  die  space-time  focusing  is  improved 
when  both  the  bandwidth  and  the  delay-spread  are  increased. 
Figs.  13  and  14  show  E{t]d(11)/7}d(&o)}  and  E{ At(R)}  as 
a  function  of  die  system  bandwidth.  The  channel  delay  spread 
is  set  to  0.8  jus  for  these  simulations. 


Fig.  13.  Simulated  energy  function  R)/?J£>(Ro)}  as  a  function  of 

range  ami  bandwidth  (trms  —  0.8/xs). 


The  simulated  dependencies  have  been  fitted  by  the  follow¬ 
ing  empirical  formulas: 


E{pd}=Pd  =  5.9t™IB0-35 


E{pm}  =  Pm  =  2-2trmsB0'52  (20) 


x  I  Jo  (2tt|R  -  Rol/A)  r  exp  (~|R  -  Ro|/3.2A) 


with 

and 


a  =  2.6exp(— trms/1-23) 

Kjy  a  1/pD 

Gal  A  =  Gx/X  =  0.41r^s01B-°02 
S7o  =  O.OSSr^fB-1  +  O.93492 


(21) 

(22) 

(23) 

(24) 

(25) 


E{a( R)}  2  0.37(1  -  |  J0(2tt|R  -  Ro\/X)f 

x  exp(-|R  —  Ro|/1.5A)}  (26) 

with  (3  —  1.3exp(-7)tMs/l-8).  Here,  the  bandwidth  B  is 
expressed  in  MHz,  and  the  delay-spread  trms  is  expressed 
in  ^s. 

The  alternative  energy  ratio  £{7?Af(R)/?7Af(Ro)}  can  be 
expressed  in  the  same  fashion  as  in  (21). 

•  In  (19),  (24)  and  (25),  exponents  for  B  and  trms  are 
similar  when  both  variables  appear  multiplied  by  each 
other.  This  observation  implies  that  space-time  focusing 
in  general  is  improved  by  increasing  the  product  B  x 
trms-  However,  that  assertion  should  be  nuanced,  since 
some  terms  in  (25)  and  (26),  as  well  as  a  and  j3  are  only 
related  to  the  channel  delay-spread. 

•  When  B  is  sufficiently  large  (>  10),  file  space-focusing 
parameters  depend  upon  B  x  trms-  This  impact  is  sig¬ 
nificant  regarding  the  asymptotic  SFG.  However,  the 
dependence  of  the  3-dB  contour  width  towards  B  x  trms 
is  much  weaker.  Both  Ga  and  Gx  are  therefore  mostly 
inversely  proportional  to  the  carrier  frequency  only.  Note 
that  for  smaller  bandwidths  (for  which  eqs.  (19)  and  (20) 
are  not  valid  anymore),  the  average  asymptotic  SFG  tends 
to  unity. 

•  The  average  delay  spread  at  the  focal  point  A r0  results 
from  the  additive  combination  of  two  terms:  one  roughly 
related  to  B  x  trms,  and  the  other  one  to  trms  only.  For 
large  B  x  trms,  the  variation  of  Ato  with  B  becomes 
weak,  so  Atq  only  depends  upon  the  channel  delay- 
spread.  This  is  clearly  visible  on  Fig.  14.  At  any  other 
point  R  outside  the  focusing  region,  the  average  delay 
spread  E{  At(R)}  is  an  oscillating  increasing  exponen¬ 
tial  function  of  the  distance  from  the  focal  point.  It 
saturates  at  long  distances. 

•  At  any  point  R,  the  time- focusing  gain  £{<r(R)}  is 
only  related  to  the  channel  delay-spread.  Furthermore, 
the  asymptotic  TFG  (i.e.  at  long  distance)  is  a  constant 
value  (~  0.37)  independent  from  B  and  trms,  at  least  for 
sufficiently  large  B  x  trms  (>  10). 


(19) 
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Fig.  14.  Simulated  RMS  delay  spread  £?{Ar(R)}  as  a  function  of  range 
and  bandwidth  (trMs  =  0.8/zs). 


3)  Impact  of  Reduced  Angular  Spread :  When  the  angle 
spread  is  decreased  from  2ix  to  smaller  angles,  the  spatial 
focusing  is  degraded  as  expected. 

»  For  Ai?  >  7r/3,  the  asymptotic  SFG  po  decreases  as  At? 
decreases,  and  is  found  to  vary  as  (At?/27r)0'6. 

•  The  3-dB  contours  Ga  and  Gx  are  not  identical  any¬ 
more,  and  Ga  —  4 Gx  for  At?  <  37r/4.  They  are 
decreased  respectively  by  At?- 1,3  and  A#"1'2  relatively 
from  their  previous  values  (i.e.  when  A#  was  equal  to 
2?r).  For  At?  >  37r/4,  the  ratio  <3a(At?)/<5a(At?  —  2n) 
continues  to  decrease  following  a  At?-13  law  while 
<5x(At?)/(5x(At?  =  2?r)  remains  constant.  Hence,  Ga  = 
Gx  for  At?  =  27 r. 

New  empirical  expressions  can  therefore  be  derived  easily 
that  take  into  account  that  die  channel  delay  spread  trms  is 
also  modified  when  reducing  the  angle  spread,  as  highlighted 
by  (15)  and  (16): 

Pd  =  5.9T&iB0'35(A^/27r)0V-°-2  (27) 

with  v  —  (A#  +  sin  Ai?)/(2tt), 

Ga/\  =  0.41r^s01B-0  02(Ai?/27r)_1'3i/0  005 

w  0.41(Ai?/2?r)-1,3  (28) 

[  0.13r^s01B-002(At?/2^)-1-2t/0005, 

Gx/\  =  \  tr/3  <  At?  <  3tr/4, 

[  0.41trj^s01B~0'<wi'0  005,  At?  >  3tr/4 

_  f  0.13(At?/2w)-1-2,  tt/3  <  At?  <  3tt/4, 

~  |o.41(At?/27r)~1-3,  At?  >  3tt/4. 

4)  Benefits  from  Multiple  Antennas:  The  impact  of  using 
multiple  antennas  at  the  transmit  point  is  analyzed  in  the 
following.  At  the  transmit  point,  the  terminal  now  consists 
of  Mt  antennas.  Each  antenna  has  a  specific  location  around 
T,  denoted  as  Tu.  Hence,  the  focused  signal  is  written  as 

s(r,  Ro)  = 

1 

— V  h{r,  Ro  ->  Ttt)  ★  h(-r,  Ro  -►  T„)  (30) 


where  the  scaling  factor  1/y/Mr  ensures  that  the  total  trans¬ 
mitted  power  is  kept  constant,  irrespective  of  the  number  of 
antennas. 

The  benefit  from  multiple  antennas  is  proportional  to  Mt 
as  far  as  the  SFG  is  concerned,  in  agreement  to  well-known 
diversity  results  [31].  Indeed,  considering  a  constant  transmit¬ 
ted  power,  the  peak  energy  at  the  focal  point  will  grow  as  Mr, 
but  will  remain  the  same  anywhere  else  outside  the  focusing 
region. 

5)  Multiple-Bounce  Model:  Single-bounce  models  are  not 
sensitive  to  the  order  of  magnitude  of  T  and  7  since  they  are 
only  responsible  for  a  proportionality  factor.  This  is  not  the 
case  for  multiple-bounce  models  in  which  the  average  loss 
per  interaction,  related  to  the  scattering  cross-section  and  the 
path-loss,  is  a  critical  parameter. 

The  model  can  be  easily  extended  to  consider  multiple 
bounce  interactions.  In  such  case,  the  transmitted  waves  are 
scattered  from  one  to  a  certain  number  of  times.  The  calcu¬ 
lation  of  the  transmission  channel  can  be  written  as  a  sum  of 
N  scattered  contributions: 

N 

H( «,Ro)  =  £Ci»)iG(w,S,  ->T)  (31) 

2—1 

where 

•  G(iu,  Si  — >  T)  is  the  spreading  function  for  a  transmis¬ 
sion  from  the  Ith  scatterer  to  the  target  point  T, 

•  771  is  the  scattering  cross-section  of  the  Ith  scatterer, 
proportional  to  the  scattering  coefficient  Ti  [33], 

•  Ci  is  die  Ith  element  of  vector  C,  which  is  the  solution 
of  a  iV-order  linear  system  given  by: 

C  =  A-1  •  Co  (32) 


with 

•  Cqi  ~  G(u;,Ro  — >  Si)  is  the  channel  between  die  target 
point  and  the  location  of  each  scatterer, 

•  A  is  the  system  matrix  accounting  for  the  transmission 
between  each  pair  of  scatterers: 


irkG(u,sk-+&i) 

\l  k  =  L 


(33) 


As  already  mentioned,  this  model  is  highly  sensitive  to  the 
average  loss  per  interaction,  i.e.  to  the  product  \rjkG(uj ,  S*  — ► 
S|)[.  On  one  hand,  if  the  latter  is  very  small,  then  the 
delay-spread  of  the  channel  tends  to  be  infinite,  which  is 
unrealistic.  On  the  other  hand,  if  the  loss  per  interaction  is 
very  large,  the  impact  of  high- order  scattered  contributions 
becomes  small,  and  the  received  field  is  very  similar  to  the 
first-order  field.  A  reasonable  assumption  is  that  contributions 
up  to  the  third  order  can  be  significant.  The  average  loss  per 
interaction  should  then  be  chosen  accordingly.  For  an  average 
loss  per  additional  reflection  of  12  dB,  the  impact  of  multiple 
reflections  on  the  focusing  is  negligible.  If  the  average  loss 
per  interaction  is  decreased  by  4  dB,  the  impact,  though  not 
negligible,  remains  small.  This  is  mostly  because  the  large 
number  of  scatterers  enables  the  time-reversed  field  to  be 
highly  focused  with  single-bounce  interactions. 
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6)  Comparison  with  Random  Waveguide  Model:  In  both 
models,  the  space-focusing  gain  increases  with  bandwidth.  In 
the  geometry-based  model,  it  is  also  shown  that  the  SFG 
increases  similarly  with  the  channel  delay-spread.  In  the 
waveguide  model,  simulation  results  illustrate  that  the  SFG 
also  increases  with  the  channel  delay-spread.  The  increase  for 
the  waveguide  is  related  to  the  source  angle  and  to  the  inverse 
of  the  waveguide  width.  Time-focusing  gain  examined  using 
the  3-dB  contour  shows  a  weak  dependence  on  bandwidth 
in  both  scenarios.  Similarly,  the  dependence  on  the  channel 
delay-spread  is  also  weak,  but  the  3-dB  contour  decreases  as 
the  angle-spread  increases.  This  is  clear  in  the  geometry-based 
model.  In  the  waveguide  model,  we  have  observed  that  the 
dependence  toward  the  waveguide  width  is  weak,  while  there 
is  a  noticeable  dependence  on  the  angle-spread. 

Some  differences  between  these  two  models  have  been 
observed.  The  most  pronounced  one  is  that  the  SFG  appears 
to  be  independent  of  the  bandwidth  in  the  case  of  the  discrete 
scattering  model,  while  a  significant  dependence  is  found  in 
the  waveguide.  Also,  the  absolute  value  of  the  SFG  is  much 
larger  in  the  waveguide  scenario.  To  better  understand  this 
difference,  we  need  to  refer  to  the  specifics  of  the  models. 
In  the  waveguide  case,  two  mechanisms  cause  delay-spread: 
large-scale  multipath  results  from  the  coherent  reflections  by 
the  waveguide  walls,  small-scale  multipath  is  caused  by  the 
non-homogeneities  of  the  refractive  index.  By  contrast,  the 
geometry-based  propagation  only  consists  of  non-coherent 
reflections,  similar  to  the  waveguide  small-scale  multipathing. 
Focusing  in  time  reversal  is  enhanced  for  bounded  domains  in 
that  it  is  more  robust  This  accounts  for  the  difference  seen  in 
these  two  models.  As  the  waveguide  width  approaches  infinity, 
it  begins  to  agree  with  an  unbounded  medium. 

7)  Space-time  focusing  in  time-reversed  transmissions:  So 
far,  we  have  assumed  that  a  single  time-reversed  pulse  is  sent 
from  the  transmitter  to  the  target  point.  However,  in  Section  V, 
we  will  deal  with  data  transmissions  in  TR  schemes.  As  a 
preliminary  study,  we  analyze  here  how  space-time  focusing 
is  affected  when  sending  a  time-reversed  data  signal.  The  latter 
consists  of  the  convolution  of  the  time-reversed  pulse  with  a 
train  of  information  bits. 

Fig.  15  shows  the  SFG  as  a  function  of  the  relative  data 
rate  (i.e.  the  data  rate  relative  to  the  pulse  bandwidth)  for  two 
different  Nyquist  roll-off  factors.  The  bandwidth  is  equal  to  40 
MHz  and  the  channel  RMS  delay-spread  to  0.8/iS.  For  these 
values  the  single-pulse  SFG  Pd  ~  13  dB.  It  is  clear  from 
Fig.  15  that  for  the  full  rate  the  space  focusing  advantage  is 
totally  destroyed.  The  same  conclusion  applies  regarding  the 
TFG,  as  depicted  in  Fig.  16.  In  other  words,  for  a  single  user, 
the  time  compression  effect  resulting  from  TR  is  significantly 
reduced  if  the  rate  is  not  seriously  decreased. 

In  fact.  Fig.  15  and  16  show  that  space  and  time  focusing 
are  slowly  restored  as  the  data  rate  decreases.  As  an  example, 
the  SFG  is  about  5  to  6  dB  for  a  rate  of  B/ 10,  while  the  TFG 
is  around  0.22  for  that  same  rate.  For  higher  roll-off  factor,  the 
degradation  is  less  than  for  the  zero  roll-off  case.  This  is  not 
surprising  since  higher  roll-off  factors  yield  lower  side  lobes. 

Therefore,  there  is  a  fundamental  trade-off  between  the  data 
rate  and  the  space-time  focusing  effect  that  can  be  achieved 


Fig.  15.  Spatial  focusing  gain  as  a  function  of  data  rate  and  roll-off  factor. 


Fig.  16.  Time  focusing  gain  as  a  function  of  data  rate  and  roll-off  factor. 

in  a  wireless  channel  for  a  given  bandwidth.  This  trade-off  is 
dictated  by  the  coherence  bandwidth  of  the  channel  transfer 
function  and  is  analyzed  in  more  detail  in  the  next  section. 

V.  Communications  over  Large  Delay  Spread 
Channels  with  Very  Simple  Receivers 
In  this  section,  we  investigate  the  performance  of  a  most 
simple  TR  communication  system.  We  assume  that  the  channel 
shortening  achieved  by  TR  is  all  that  is  done  to  combat 
remaining  ISI  in  the  channel  and  investigate  an  uncoded 
MISO  system  that  employs  a  simple  one-tap  receiver  without 
any  further  equalization.  We  study  the  performance  of  this 
TR  system  in  the  limit  of  very  low  and  of  infinite  delay 
spread.  We  compare  this  performance  to  that  of  an  optimized 
preequalizer,  and  examine  the  impact  of  the  number  Mt  of 
transmit  antennas,  and  of  a  strategy  called  rate  back-off  that 
can  combat  remaining  ISI  in  the  channel  in  a  very  simple  way. 
The  goal  of  this  study  is  two- fold.  We  obtain  a  lower  bound  on 
the  performance  of  a  TR  communication  system  in  channels 
with  very  high  delay  spread,  since  any  other  system  that  does 
employ  coding  or  equalization  at  the  receiver  will  perform 
better.  Note  that  the  upper  performance  bound  was  presented 
in  Section  III,  since  the  computation  of  the  system’s  rate  in 
that  section  assumed  both  infinitely  complex  transmitter  and 
receiver  and  perfect  coding.  Secondly,  the  proposed  setup  can 
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Fig.  17.  A  block  diagram  of  a  MISO  system  with  a  preequalizer  at  the 
transmitter  and  a  single  tap  receiver. 


be  of  interest  by  itself  for  systems  that  operate  with  very  large 
banwidths,  not  too  high  data  rates,  and  that  rely  on  extremely 
simple  receiver  architectures. 

Fig.  17  depicts  a  block  diagram  of  the  system.  A  symbol 
x  is  upsampled  by  a  factor  of  D ,  put  on  the  Mt  branches 
to  the  transmit  antennas,  filtered  with  a  function  gz  which  can 
possibly  be  a  TR  prefilter,  and  sent  over  the  channels  hx.  After 
receive  noise  is  added,  the  signal  can  be  enhanced  by  a  gain 
/?,  and  it  is  downsampled  before  detection. 

For  all  plots  shown  in  this  section,  we  use  a  Rayleigh  fading 
channel  with  L— 100  i.i.d.  taps. 

A.  Optimal  Prefiltering 

First  we  derive  the  optimal  prefilter  for  the  system  in 
Fig.  17.  We  describe  the  input-output  relation  of  the  system 
in  matrix  form  as 

y[k]  =  /?xfcHg  +  (3n[k ]  (34) 

where  xfc  =  [xix2  . .  .^22,-1]  is  a  vector  of  input  symbols, 
y[k\  the  output  symbol  and  n[k]  the  noise.  The  discrete  time 
channel  matrix  H  is  block  toeplitz;  its  block-columns  are 
comprised  of  the  shifted  channel  vectors  [hi  [1], . . . ,  hi[L]],  and 
in  each  block  columns  that  belong  to  adjacent  antennas  are 
next  to  each  other.  The  vector  g  is  given  by 

g  =  [g[0]r  ...  slL-lff  (35) 

and 

m  =  [9i[l\  ...  9mt[1]\T  /  =  (36) 

The  product  Hg  is  die  convolution  of  the  prefilter  and  the 
channel;  T  denotes  die  transpose. 

We  give  here  the  optimal  g  and  j3  for  the  case  of  D  —  1 
(no  upsampling).  In  order  to  simplify  receiver  processing, 
the  purpose  of  the  set  of  transmit  filters  g  is  to  preequalize 
the  signal.  This  means  that  at  the  receiver,  as  much  power 
as  possible  is  focused  in  the  central  peak,  and  that  ISI  is 
minimized,  subject  to  the  power  constraint  ||  g  ||=  1  where 
||  .  ||  denotes  the  norm.  Such  a  problem  can  be  stated  as  an 
optimization  problem,  whose  solution  is 

*-KH"H+Ki)''H"'a  <37> 

/J2  =  e^H  (hhH  +  j^-lj  HHeA.  (38) 

A  is  the  delay  of  the  equalizer,  and  the  unit  vector  with 
1  in  entry  A,  and  Ptx  —  £-[xfcxfc]*  Eq.  (37)  defines  a  filter 


Fig.  18.  The  focused  CIR  for  the  TR-  and  the  optimal  prefi  Iter. 


with  an  MMSE-type  solution  with  a  finite  number  of  taps.  We 
notice  that  the  optimized  value  for  f3  bascially  ensures  that  the 
power  constraint  of  the  filter  is  valid.  In  the  transmit-receiver 
chain  of  Fig.  17,  0  has  no  effect  on  the  SNR  at  the  receiver;  it 
ensures,  however,  that  the  transmitter  can  take  into  account  the 
SNR  to  optimize  the  preequalizer.  We  note  that  even  though 
there  exists  a  striking  similarity  to  an  MMSE  receiver,  we 
could  not  establish  the  convexity  of  the  optimization  problem. 
Hence,  our  solution  (37)  may  actually  only  be  a  local  optimum. 
The  similarity  of  the  solution  to  the  MMSE  receive  equalizer 
suggests  that  the  found  optimum  shows  good  performance;  but 
better  solutions  may  exist. 

The  TR-prefilter  is  a  matched  prefilter;  it  follows  from  the 
optimal  filter  (37)  for  low  SNR, 


gTR 


lim 

SNR-+0 


g  = 


h* 


s.t.  II  h  11=1, 


(39) 


or 

with  7  =  (Et=i  Ef'Jo1  IftiWI2)'1.  i.e.,  the  filter  equals  the 
time-reversed  CIR  of  the  channel,  and  the  receiver  receives 
the  autocorrelation  function  of  the  channel.  The  power  in  its 
central  peak  is  maximized,  but  the  ISI  suppression  is  worse. 

A  picture  of  the  CIRs  that  result  from  precoding  with  the 
optimal  and  the  TR  strategy  is  shown  in  Fig.  18.  The  SNR 
is  10  dB.  The  optimal  prefilter  suppresses  the  ISI  better  than 
the  TR  prefilter;  close  inspection  reveals  that  the  power  in 
the  peak  of  the  CIR  is  slightly  higher  in  case  of  TR.  The 
improved  ISI  suppression  in  case  of  die  optimal  prefilter  comes 
at  computational  cost  The  optimal  transmit  prefilter  requires 
a  matrix  inversion  which  can  be  difficult  if  the  number  of 
channel  taps  is  too  high  or  the  channel  is  not  sufficiently 
stationary. 

B.  Rate  Back-off 

Fig.  1 8  suggests  that  the  application  of  TR  or  even  of  an 
optimized  prefilter  with  a  single  antenna  leaves  considerable 
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residual  I  SI  in  the  channel.  One  method  to  combat  ISI  is  to 
transmit  at  lower  rate;  if  the  transmission  rate  is  by  a  factor 
of  D  lower  than  the  bandwidth  of  the  channel,  the  ISI  is 
on  average  also  reduced  by  a  factor  of  D ,  but  the  useful 
power,  concentrated  in  the  peak  of  the  CIR,  remains  constant. 
Such  a  strategy  is  not  favorable  if  high  spectral  efficiency 
is  required.  But  for  systems  which  operate  at  very  large 
bandwidths,  such  as  ultra-wideband  systems,  it  can  be  very 
attractive  to  trade  abundantly  available  transmission  bandwidth 
against  complexity. 

Rate  back-off  by  a  factor  of  D  is  equivalent  to  upsampling 
the  signal  before  transmission,  and  downsampling  it  by  D 
before  detection.  The  input-output  relation  of  the  pre-equalized 
channel  with  rate  back-off  factor  D,  D  =  1, . . . ,  L>  can  be 
derived  as 

Mr 

y^{k]=^J2^hi*9i)mx[k-l}  +  Pn{k}  (40) 

1=1  l 


where  p  =  . 

We  furthermore  define  the  SNR  pmfb  that  would  be 
obtained  if  there  was  no  precoding  but  the  receiver  had  a 
matched  filter  and  used  full  channel  knowledge  as 

Pmfb  =  ~~ - - •  (4o) 

2)  Application  to  TR  Systems:  The  simplicity  of  the  TR 
prefilter  allows  to  compute  the  ratio  k  of  the  useful  power 
against  the  ISI  for  an  exponentially  decaying  channel  with 
delay  spread  err.  One  obtains  in  the  limit  of  large  L/or  [29] 

kTr  =  ^1  —  exp  — 

[Mt(  1  +  exp  -  i)  +  MT(1  -  exp  -  i)j 
2MTexp -£•  (l  —  exp-^) 

(47) 


where  x[k]  is  the  input  symbol  and  y^[k]  the  output  symbol 
received  with  rate  back-off.  We  see  that  the  rate  back-off  yields 
an  equivalent  channel  which  is  sampled  only  at  each  Dth 
instant. 


C.  Focusing  Ratio  and  Effective  SNR 

1)  Definitions:  We  introduce  some  performance  measures 
that  we  will  later  use  for  system  characterization.  A  very 
intuitive  performance  measure  is  the  focusing  ratio  k  which 
we  define  as  .  . 

EM  (4i) 

*  -  E[QV>  1]  (  ’ 


where  E[PQ]  is  the  expected  power  in  the  central  peak  of 
the  preequalized  CIR,  and  E[Q^\  is  that  of  the  ISI  with  rate 
back-off  D,  i.e.,  the  power  in  the  tails  of  the  preequalized  CIR. 
The  better  the  focusing,  the  larger  is  k.  Unlike  the  measure 
introduced  in  (9),  the  measure  here  takes  into  account  the 
entire,  temporally  spread  power  of  the  received  signal,  and 
not  just  the  ratio  of  the  energies  of  the  maximum  and  of  one 
distant  tap. 

We  can  express  the  power  of  the  received,  possibly  down- 
sampled,  signal  as 

£[|y(DW]  -  PTxE[P0\+PrxE{Q[n]}  +  a2  (42) 


where  we  used  E[x[l]x[m]]  =  Prx<5[/— m],  with  Ptx  denoting 
the  transmitter  power,  and  E[n[l]n*[m]]  =  o2S[l  -  m }  with 
cr2  the  noise  power  and  6  [l  —  m]  the  Kronecker  Delta.  The 
constants  Po  and  are  given  as  [29] 

j  Mr  L—l  I2 


P 


EE  hi[l}gi[-i] 


(43) 


qm  = 


EEE  hi[l}gi[m}6[Dk 


i=  i  &o  k^o 


—  m  — 1\ 


2 

1  .(44) 


We  define  the  effective  SNR  pej  /  as 


PtxE[P0) 

Peff,D  PtxE[QM]  +  a2 


E\Pq\ 

mm + ^ 


(45) 


Using  (47)  together  with  (45)  and  (46)  we  arrive  at  an  analyt¬ 
ical  expression  for  the  effective  SNR  of  a  TR  communication 
system  with  Mr  transmit  antennas  and  a  rate  back-off  of  P, 

peff.D  =  PMFB  <48) 

'  «Ti? 

We  can  derive  from  this  the  following  limiting  cases: 


Peff,D 


Pmfb  for  or  — ►  0 

PmfbDMt 

- -  -  -  ror  oT 

Pmfb  +  DMt 


(49) 
oo  (50) 


Eq.  (50)  can  be  split  into  a  low-SNR  and  a  high-SNR  regime 
as 


Pefj,o  =  Pmfb  for  pmfb  low,  oT  large  (51) 
Peff,D  =  DMt  for  pmfb,  crT  large.  (52) 

Interpretation  of  these  equations  yields  the  following  conclu¬ 
sions: 

•  Eq.  (49)  tells  that  a  TR  system  is  lossless  if  the  delay 
spread  of  the  channel  is  very  small.  Preequalization  of  the 
channel  does  hardly  change  the  performance  if  the  delay 
spread  is  very  small.  We  recall  from  Fig.  2,  however, 
that  the  system  in  this  case  is  subject  to  fading  which 
will  reduce  the  rate. 

•  Eq.  (51)  shows  that  TR  is  optimal  at  low  SNR;  this 
behavior  is  expected  for  a  transmit  matched  filter  which 
maximizes  the  power  in  the  peak  of  the  signal. 

•  Eq.  (52)  demonstrates  that  in  an  environment  with  a  large 
delay  spread,  the  performance  of  a  TR  system  at  high 
input  SNR  saturates  at  a  value  independent  from  the  delay 
spread  of  the  channel.  The  performance  increases  with 
3  dB  per  added  transmit  antenna,  and  with  another  3  dB 
for  each  rate  back-off  by  a  factor  2  from  the  maximum 
available  bandwidth. 

The  last  property  is  remarkable,  because  it  holds  for  an 
uncoded  system  with  no  receive  equalizer,  a  fairly  simple 
transmit  equalization  strategy,  and  arbitrary  delay  spread,  i.e., 
arbitrary  bandwidth.  In  a  1  GHz  time-invariant  channel,  a  sys¬ 
tem  with  a  powerful  two-antenna  transmitter  can  theoretically 
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Fig.  19.  The  effective  SNR  increases  with  the  number  of  transmit  antennas.  Fig.  20.  Rate  back-off  improves  the  effective  SNR. 


achieve  a  transmission  rate  of  125  Mbps  at  12  dB  receive 
SNR  with  a  single  one-tap  receiver  and  independent  of  the 
delay  spread  of  the  channel. 


D.  Performance  Evaluation 

We  study  the  performance  of  the  optimal  and  the  TR 
prefilters  in  dependence  on  the  number  of  transmit  antennas 
and  the  rate  back-off,  and  examine  furthermore  the  variance 
of  the  received  power  as  a  function  of  bandwidth. 

In  Fig.  19  we  plot  the  effective  SNR  of  the  two  filters 
as  a  function  of  pmfb  ranging  from  -10  to  20  dB;  in  this 
picture,  an  optimal  system  would  be  characterized  by  a  line 
with  slope  1.  The  effective  SNR  is  computed  using  100  CIRs. 
The  curves  for  die  MMSE  filter  are  dashed,  those  for  the  TR 
filter  continuous.  The  three  different  curves  for  each  filter  are 
obtained  for  Mt  =  1,2,  and  4.  The  results  for  the  TR  filter 
match  the  prediction  of  (47).  For  large  pmfb >  the  effective 
SNR  becomes  constant.  It  is  0  dB  for  Mt  =  1,  and  increases 
with  3  dB  for  each  additional  transmit  atenna.  At  low  SNR, 
the  curves  increase  linearly  with  pmfb *  In  this  range,  the 
TR  filter  performance  equals  that  of  the  MMSE  filter.  The 
effective  SNR  for  the  SISO  system  is  quite  bad;  however,  the 
shown  performance  is  achieved  with  an  extremely  simple,  one- 
tap  receiver,  and  in  a  channel  with  very  long  delay  spread. 
The  optimized  prefilter  shows  much  better  performance.  For 
Mt  =  1,  the  effective  SNR  saturates  at  a  value  which 
appears  to  depend  on  the  delay  spread  of  the  channel  and  its 
correlation  between  taps.  For  more  than  one  transmit  antenna, 
the  optimized  prefilter  can  completely  eliminate  ISI,  so  that 
peff  grows  without  bound.  The  slope  of  the  curve  depends  on 
the  number  of  transmit  antennas.  We  believe  it  is  interesting 
to  study  analytical  expressions  for  the  behavior  of  pejf  for 
the  optimal  prefilter  in  the  limit  of  very  large  bandwidth. 

The  system  can  operate  with  rate  back-off  in  order  to 
increase  the  performance  at  the  cost  of  transmission  rate.  The 
effective  SNR  for  a  single-antenna  system  with  rate  back-off 
of  D  =  1 , 2, 5,  and  25  is  displayed  in  Fig.  20.  It  increases  with 
decreased  rate.  Also,  the  saturation  of  the  effective  SNR  sets  in 


much  later;  the  (almost)  linear  increase  of  pej  f  reaches  for  the 
TR  system  with  D  =  5  up  to  pmfb  =  10  dB,  and  for  one  with 
D  —  25  up  to  15  dB.  In  this  linear  range,  the  performance 
of  the  TR  prefilter  is,  in  terms  of  the  effective  SNR,  partly 
higher  than  that  of  the  MMSE  prefilter,  since  the  latter  is 
optimized  for  ISI  suppression  at  D  =  1  and  has  less  power  in 
its  central  peak.  If  the  MMSE  prefilter  had  been  optimized  for 
transmission  with  rate  back-off,  it  would  be  at  least  as  good 
as  the  TR  prefilter.  The  design  of  such  a  prefilter  could  be 
interesting  for  preequalized  wideband  data  transmission.  When 
spectral  efficiency  has  not  highest  priority,  transmission  rate 
back-off  (while  keeping  the  transmission  bandwidth  constant) 
is  very  efficient  in  order  to  improve  system  performance. 

We  characterize  the  fading  of  the  channel,  which  we  com¬ 
pute  as  the  normalized  variance  of  the  instantaneous  effective 
SNR, 


aSNR  z 


■-  Var- 


PtxPq  ,  pr  PtxPq  -|2 

PTxQM  +  <72/  [PTxQM+(r2l‘ 


(53) 


We  compute  the  variance  for  both  filters  for  p  —  10  dB  and 
Mt  —  1  in  dependence  on  the  number  of  taps  of  the  channel. 
If  the  channel  has  only  a  few  taps  the  peak  of  the  preequalized 
CIR  fades  with  higher  probability  than  at  high  bandwidth  with 
a  large  number  of  taps.  The  result  is  shown  in  Fig.  21  for 
the  number  of  taps  in  the  range  between  about  1  and  100. 
The  shown  curves  are  based  on  an  average  over  100  different 
realizations  of  the  channel.  We  see  that  the  variance  decreases 
quickly  with  the  number  of  taps;  at  only  20  taps  it  is  already 
very  low. 


E.  Bit-Error-Rates  of  the  System 

We  plot  the  bit-error-rates  (BERs)  of  the  system,  for  rate 
back-offs  of  D  =  1,2,5,  and  25,  i.e.  the  same  as  in  Fig.  20. 
In  Fig.  22  we  show  the  BER  of  the  TR  (blue,  cont.)  and 
the  MMSE  (black,  dashed)  prefilter  against  SNR.  Simulations 
were  obtained  from  sending  5000  bits  over  5000  different 
channels  with  L  —  100.  For  D  =  1  and  also  D  —  2,  BER  of 
the  TR  prefilter  floors  quickly.  For  large  rate  back-off,  BER  is 
due  to  the  increased  effective  SNR  in  the  same  range  as  that 
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Fig.  21.  The  variance  of  the  received  signal  decreases  with  increasing 
bandwidth. 


Fig.  22.  The  BER  of  the  TR  and  the  MMSE  prefi  Iter. 

of  the  MMSE  prefilter  for  lower  rate  back-off.  The  TR  system 
performs  well  for  D  >  5;  the  MMSE  prefilter  performs  well 
for  D  >  2.  For  large  D,  both  systems  perform  in  general 
similarly  well,  since  the  MMSE  filter  was  not  optimizted  for 
rate  back-off.  A  worthwhile  study  would  be  the  performance 
of  a  system  equipped  with  a  short  MMSE  equalizer  at  the 
receiver.  The  performance  would  potentially  be  much  better 
and  the  rate  back-off  could  be  reduced  even  though  only  a 
simple  receive  equalizer  was  used. 

VI.  TR  Channel  Sounding 

The  value  of  TR  for  practical  systems  can  only  be  assessed 
through  experiments.  In  this  section  we  describe  how  a  channel 
sounder  can  be  designed  that  can  perform  real-world  TR 
experiments.  This  design  does  not  require  new  technological 
concepts;  we  rather  want  to  show  how  a  conventional  wide¬ 
band  MIMO  sounder  can  be  converted  into  a  sounder  for  TR 
experiments.  A  more  detailed  version  of  the  following  can  be 
found  in  [35]. 


A.  Architecture 

In  TR  experiments,  one  wants  to  demonstrate  temporal  and 
spatial  focusing  achieved  in  a  wide  band  channel.  In  order  to 
demonstrate  temporal  focusing,  high  resolution  of  the  received 
CIR  is  required.  To  reveal  spatial  focusing,  more  than  one 
receive  antenna  is  necessary;  one  receive  antenna  is  the  one 
that  receives  the  time  compressed  signal.  At  any  other  antenna, 
a  noise-like  signal  will  be  visible.  A  wideband  SIMO  system 
is  hence  the  minimum  requirement  to  demonstrate  TR.  The 
key  difference  between  between  MEMO  channel  sounding  and 
a  TR  experiment  is  that  die  latter  requires  CSI.  The  CSI 
must  be  updated  in  intervals  the  length  of  which  depend  on 
the  Doppler  spread  of  the  channel.  Since  even  comparably 
stationary  channels  as  those  for  fixed  wireless  access  can 
have  a  Doppler  spread  of  more  than  1  Hz  [36],  the  CSI 
must  be  updated  much  faster  than  once  per  second  and  slow 
methods,  like  a  frequency  domain  sweep  measurement,  are  not 
useful.  We  will  consider  time-  or  frequency  domain  correlation 
methods,  since  they  are  both  fast  and  can  be  used  for  wideband 
channel  sounding. 

In  common  correlation-based  wideband  channel  sounding 
experiments,  a  transmitter  transmits  known  sequences  to  a 
receiver.  Transmitter  and  receiver  are  synchronized,  and  the 
receiver  uses  the  received  sequences  in  order  to  estimate 
the  channel.  In  a  TR  experiment,  the  transmitter  knows  the 
channel,  and  convolves  the  transmit  sequence  with  the  time 
reversed  CIR  before  transmission.  The  receiver  uses  again 
the  received  sequences  to  estimate  the  time-reversed  channel, 
which  will  now  exhibit  spatial  and  temporal  focusing.  In 
order  to  convert  a  MIMO  channel  sounder  into  a  TR  channel 
sounder,  CSI  must  be  made  available  to  the  transmitter,  and 
the  transmitter  must  be  able  to  convolve  the  CIR  with  the 
transmit  sequence.  Furthermore,  the  sounder  must  still  be  able 
to  estimate  the  channel  before  CSI  is  fed  into  the  transmitter. 

A  high-level  sketch  of  a  possible  design  is  shown  in  Fig.  23. 
A  MIMO  channel  sounder,  which  consists  of  an  RF  part  for 
transmitter  and  receiver  (RFTx  and  RFRx),  a  signal  generator 
at  the  transmitter  (SGTx)  and  a  data  storage  unit  at  the 
recevier  (DSRx)  is  equipped  with  a  TR-signal  processing 
unit  (TRSPU).  The  main  task  of  the  TRSPU  is  to  drive  the 
estimiation  of  the  CSI,  and  to  transfer  the  CSI  via  the  feedback 
loop  to  the  transmitter,  where  the  CIRs  are  convolved  with 
the  transmit  sequence  and  transmitted  over  the  channel.  The 
details  of  the  TRSPU  and  the  feedback  loop  are  explained  in 
the  next  subsection. 

B.  TRSPU  and  Feedback  Loop 

A  TR  channel  sounder  operates  in  three  steps,  which  are 
controlled  by  the  TRSPU  and  shown  in  Fig.  24.  In  a  first  step, 
the  channel  is  estimated  between  a  given  transmit  antenna, 
say  Txo,  and  a  given  receive  antenna,  Rxo.  This  can  be  done 
by  sending  the  known  sequence  T(oj),  here  expressed  in  fre¬ 
quency  domain,  rii  times  over  the  channel.  Once  the  channel  is 
estimated  with  some  predefined  accuracy,  the  estimated  CSI 
#0,0  is  sent  via  the  feedback  loop  to  the  transmitter.  In  a 
third  step,  the  known  transmit  sequences  T(u)  are  multiplied 
with  the  estimated  transfer  function  and  transmitted  into  the 
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Fig.  23.  A  TR  channel  sounder  can  be  designed  using  a  MIMO  channel 
sounder  with  a  feedback  loop  for  CSI. 


Fig.  24.  The  operational  steps  of  a  TR  channel  sounder. 

channel.  In  this  third  step,  the  sounder  switches  between  all 
its  receive  antennas.  At  Rx 0,  it  will  receive  a  time  reversed 
CIR;  at  all  other  antennas,  it  receives  a  noise-like  signal.  This 
time-reversed  sequence  can  be  sent  n 2  times,  until  the  CSI 
is  outdated  and  the  sounder  has  to  start  again  with  channel 
estimation. 

Particular  care  is  required  for  the  choice  of  the  feedback 
loop.  The  transferral  time  of  die  CSI  must  be  much  shorter 
than  the  coherence  time  of  the  channel.  The  data  rate  of  the 
link  becomes  hence  a  critical  parameter.  A  single  complex  CIR 
with  about  1000  taps  and  8  bit  resolution  has  a  size  of  about 
16  kbit.  Assuming  a  channel  with  Doppler  frequency  of  1  Hz, 
a  realistic  transmission  time  should  not  be  longer  than  0.2  s; 
consequently,  the  feedback  loop  needs  a  data  rate  of  at  least 
80  kbit/s  per  CIR  that  is  transferred. 

The  following  options  can  be  considered 

1)  Microwave  Link:  A  microwave  link  is  advantageous  in 
terms  of  data  rate,  which  is  in  the  order  of  Mbit/s,  and  in  terms 
of  flexibility.  Microwave  links  can  be  established  whereever 
the  experiment  is  to  be  performed.  However,  TR  experiments 
are  likely  to  be  conducted  in  environments  with  large  delay 
spreads.  The  TR  sounder  will  operate  at  maximum  allowed 
EIRP.  This  limit  holds  for  the  microwave  link,  as  well.  In 


areas  where  the  SNR  of  the  TR  experiment  is  low,  that  of 
the  microwave  feedback  link  is  low,  too.  The  link  quality  can 
be  improved  by  using  directional  antennas  mounted  on  high 
poles,  but  this  requires  difficult  adjustments  and  large  efforts 
in  the  preparation  of  the  experiment. 

2)  Cellular  Systems:  Currently  available  cellular  systems 
such  as  GPRS  offer  great  flexibility  without  die  need  for 
heavy  experimental  equipment.  Their  main  drawback  is  their 
comparably  low  data  rate  (order  of  100  kbit/s  for  GPRS)  and, 
in  particular,  the  unreliable  transmission  delay.  Since  the  data 
rate  offered  does  not  have  too  much  overhead  compared  to  the 
data  rate  required,  the  application  of  current  cellular  standards 
for  the  feedback  loop  appears  risky  in  particular  if  the  Doppler 
spread  of  the  channel  may  at  times  exceed  1  Hz. 

3)  Ethernet:  The  ethemet  offers  tremendous  data  rates;  its 
main  drawback  is  the  inflexibility  of  its  access.  However,  most 
TR  experiments  will  be  conducted  in  urban  areas,  where  access 
to  the  ethemet  can  be  expected  to  be  in  range.  A  wireless 
network  system  such  as  one  based  on  the  IEEE  802.11b 
standard  still  offers  data  rates  in  the  order  of  Mbits  and  can  be 
employed  to  bridge  the  gap  between  the  experimental  device 
to  the  next  ethemet  connector. 

The  final  decision  about  the  best  link  depends  on  the 
availability  of  particular  services,  on  local  network  traffic,  and 
the  Doppler  spread  of  the  channel. 

C.  Further  Critical  Specifications  and  Sources  of  Error 

1)  Range  and  Power:  TR  works  best  in  environments  with 
large  delay  spreads,  i.e.,  large  distances.  Consequently,  the 
demands  on  the  link  budget  are  high.  Critical  parameters 
which  determine  range  are  here  pathloss  exponent,  maximum 
peak  and  average  transmit  power  determined  by  hardware  and 
regulations,  bandwidth  and  carrier  frequency,  antenna  gains, 
correlation  gains,  the  sounder’s  sensitivity  and  target  SNR. 

2)  Phase  Noise:  TR  requires  exact  phase  information  at 
the  transmitter.  Inaccuracies  in  the  phases  occur  since  a)  the 
channel  is  not  perfectly  time  variant,  b)  the  LO  signals  derived 
by  PLLs  from  the  (usually)  10  MHz  reference  have  phase 
noise.  The  time  variance  can  be  overvome  if  the  data  rate  of 
the  feedback  loop  is  large  enough. 

Phase  noise  destroys  the  coherent  superposition  of  signals 
transmitted  from  multiple  antennas.  Since  it  enters  both  at  the 
transmitter  and  the  receiver,  it  will  reduce  the  received  power 
by  a  factor  of  cos(a<£i  +  a <£2),  where  a 4>u  i  —  1,2  are 
random  variables  that  describe  the  phase  fluctuations  of  the 
transmitter’s  and  the  receiver’s  oscillators.  For  phase  errors 
less  than  15°,  this  yields  a  performance  loss  of  about  10  %  in 
SNR,  which  is  tolerable. 

VII.  Conclusions  and  Outlook  on  Future  Work 

In  this  paper  we  introduced  the  concept  and  goals  of  time 
reversal  (TR)  techniques  in  the  context  of  broadband  wireless 
communications.  TR  by  itself  is  a  very  simple  preequalization 
scheme  which  moves  the  system’s  matched  filter  into  the 
radio  channel.  When  applied  to  wideband  channels  with  very 
large  delay  spreads,  this  preequalizer  very  effectively  captures 
the  particular  advantages  of  this  channel,  i.e.  temporal  and 
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spatial  focusing  as  well  as  channel  hardening.  TR  allows  to 
shift  complexity  from  transmitter  to  receiver  and,  thanks  to 
its  comparably  low  computational  complexity,  is  likely  to  be 
a  feasible  scheme  even  if  original  channel  impulse  responses 
are  very  long. 

We  showed  that  the  average  rate  of  a  TR  scheme  at  low 
SNR  is  very  close  to  the  ergodic  capacity  of  the  frequency 
selective  channel.  In  fading  channels,  the  outage  capacity  can 
be  improved  if  one  transmits  at  higher  bandwidth  and  exploits 
the  channel  hardening  provided  by  TR. 

Temporal  and  spatial  focusing  were  investigated  with  three 
examples  based  on  mesaured  indoor  UWB  data,  and  on  models 
of  discrete  and  of  continuous  scattering.  We  illustrated  how 
several  proposed  characterization  metrics  vary  as  a  function  of 
system  and  channel  parameters,  such  as  the  bandwidth,  delay 
and  angle  spreads,  and  number  of  antennas. 

When  addressing  communication  theoretic  issues,  only  the 
most  simple  communication  scheme  with  a  single  tap  receiver 
has  been  considered.  We  saw  that  such  a  system  can  work  in 
this  harsh  channel  when  the  transmission  bandwidth  is  sacri- 
fied  for  temporal/spatial  focusing.  More  optimal  transmission 
methods  can  achieve  much  better  performance  at  the  cost  of 
more  complexity  at  the  transmitter. 

Finally,  we  outlined  how  the  TR  channel  can  be  measured 
with  a  system  which  is  basically  an  upgraded  commercially 
available  MIMO  channel  sounder. 

This  overview  touches  many  issues  about  transmission 
over  channels  with  large  delay  spread  for  which  TR  is  one 
potentially  useful  transmission  scheme.  Many  questions  have 
not  been  addressed.  TR  is  optimal  at  very  low  SNR  for  very 
simple  receivers.  A  natural  question  is  to  ask  what  happens  if 
the  receiver  has  more  than  just  one  tap,  but  still  far  not  enough 
to  equalize  the  channel  on  its  own.  It  is  not  obvious  that  this 
problem  results  in  a  MMSE  receiver,  which  could  be  used  in 
combination  with  the  TR  precoder. 

No  space  in  this  paper  has  been  devoted  to  the  analysis  of  a 
communication  system  that  makes  use  of  the  spatial  focusing 
properties  of  the  proposed  preequalization  schemes.  This  topic 
is  inherently  connected  to  multiuser  communications.  The  rate 
region  of  such  a  system  requires  the  study  of  fading,  frequency 
selective  broadcast  channels  in  the  limit  of  large  bandwidths. 
The  derivation  of  optimal  transmit  and  receive  schemes  would 
continue  the  work  presented  in  this  paper,  or  related  work 
about  joint  transmit/receive  optimization  for  multiuser  commu¬ 
nications  but  in  the  case  of  large  delay  spread.  For  TR,  there 
appears  to  be  a  natural  relation  to  the  interference  channel.  In 
TR,  the  channel  impulse  response  is  temporally  and  spatially 
compressed  solely  with  the  knowledge  of  the  channel  of  one 
individual  link;  no  knowledge  about  the  link  to  any  other  user 
or  network  is  required.  Hence,  another  potential  application 
of  TR  is  wideband  networking  with  strong  access  points, 
architecturally  very  simple  users,  and  a  very  low  degree  of 
cooperation. 

A  very  critical  aspect  that  has  only  briefly  been  discussed 
in  Section  VI  is  that  of  channel  estimation.  We  believe  that 
TR  or  any  other  preequalized  system  can  be  applied  in  TDD 
systems  where  the  channel  state  information  is  estimated  at  the 
base  station  in  the  uplink,  and  then  used  in  the  downlink.  That 


way,  only  the  powerful  base  station  needs  to  fully  estimate  the 
channel,  and  the  handsets  only  estimate  as  many  coefficients 
as  the  short  length  of  their  receive  equalizers  requires.  In  the 
uplink,  no  preequalization  is  performed,  and  the  matched  or 
any  other  receive  filter  is  not  in  the  channel,  but  in  the  receiver 
itself. 

We  believe  that  die  study  of  communications  over  large 
bandwith  /  high  delay  spread  channels  is  very  fruitful.  These 
channels  offer  potentially  large  rate,  high  diversity,  and  spatial 
and  temporal  focusing.  TR  systems  that  operate  here  provide 
a  low-complexity  solution  for  data  transmission;  their  inves¬ 
tigation  may  be  a  suitable  starting  point  for  more  research  in 
this  area. 

Acknowledgement 

The  authors  thank  Intel  Corp.  /  Santa  Clara  (CA,  USA) 
for  providing  the  measurement  data  used  in  the  time  reversal 
project.  Claude  Oestges  is  also  pleased  to  acknowledge  the 
financial  support  of  the  Belgian  NSF. 

References 

[1]  A.  Derode,  P.  Roux,  and  M.  Fink,  ‘Robust  acoustic  time  reversal  with 
high-order  multiple  scattering”  Phys.  Rev.  Letters ,  vo!.  75,  pp.  4206- 
4209,  1995. 

[2]  M.  Fink,  ‘Time  reversed  acoustics,*  Physics  Today ,  pp.  34-40,  1997. 

[3]  M.  Fink,  ‘Time-reversed  acoustics,”  Scientifi  c  American ,  pp.  91-97,  1999. 

[4]  L.  Borcea,  G.  Papanicolaou,  and  C.  Tsogka,  ‘Theory  and  applications  of 
time  reversal  and  interferometric  imaging,”  to  appear  in  Inverse  Problems, 
2004. 

[5]  D.  R.  Dowling  and  D.  R.  Jackson,  ‘Phase  conjugation  in  underwater 
acoustics,”  J.  A  coast.  Soc.  Am .,  vol.  89,  pp.  171-181,  1990. 

[6]  D.  K  Dowling  and  D.  R.  Jackson,  ‘Narrow-band  performance  of  phase- 
conjugate  arrays  in  dynamic  random  media,”  J.  Acoust.  Soc.  Am.,  vol.  91, 
pp.  3257-3277,  1992. 

[7]  D.  R.  Dowling,  ‘Acoustic  pulse  compression  using  passive  phase- 
conjugate  processing,”/.  Acoust.  Soc.  Am.,  vol.  95,  pp.  1450-1458,  1994. 

[8]  W.  A.  Kuperman,  W.  S.  Hodgkiss,  H.  C.  Song,  T.  Akal,  C.  Feria,  and  D. 
R.  Jackson,  ‘Phase  conjugation  in  die  ocean:  experimental  demonstration 
of  an  acoustic  time-reversal  mirror”  J.  Acoust.  Soc.  Am.,  vol.  103,  pp. 
25-40,  1998. 

[9]  W.  S.  Hodgkiss,  H.  C.  Song,  W.  A.  Kuperman,  T.  Akal,  C.  Feria  and  D.  R. 
Jackson,  ‘A  long-range  and  variable  focus  phase-conjugation  experiment 
in  shallow  water,”/  Acoust.  Soc.  Am.,  vol.  105,  pp.  1597-1604,  1999. 

[10]  S.  Kim,  G.  F.  Edelmaim,  W.  A.  Kuperman,  W.  S.  Hodgkiss,  and  H. 

G.  Song,  ‘Spatial  resolution  of  time-reversal  arrays  in  shallow  water”  J. 
Acoust.  Soc.  Am.,  vol.  110,  pp.  820-829,  2001. 

[11]  D.  Rouseff,  D.  R.  Jackson,  W.  L.  J.  Fox,  C.  D.  Jones,  J.  A.  Ritcey  and 
D.  R.  Dowling,  ‘Underwater  Acoustic  Communication  by  Passive-Phase 
Conjugation:  Theory  and  Experimental  Results,”  IEEE  J.  Ocean.  Eng., 
vol.  26,  pp.  821-831,  2001. 

[12]  M.  G.  Heinemaim,  A.  Larazza  and  K.  B.  Smith,  ‘Acoustic  communica¬ 
tions  in  an  enclosure  using  single-channel  time-reversal  acoustics  ”  Appl. 
Phys.  Lett.,  vol.  80,  pp.  694-696,  2002. 

[13]  G.  F.  Edehnann,  T.  Akal,  W.  S.  Hodgkiss,  S.  Kim,  W.  A.  Kuperman  and 

H.  C.  Song,  ‘An  initial  demonstration  of  underwater  acoustic  communi¬ 
cations  using  time  reversal  ”  IEEE  J.  Ocean.  Eng.,  vol.  27,  pp.  602-609, 
2002. 

[14]  P.  Roux,  B.  Roman  and  M.  Fink,  ‘Time-reversal  in  an  ultrasonic 
waveguide,”^/?/?/.  Phys.  Lett.,  vol.  40,  pp.  1811-1813,  1997. 

[15]  P.  Roux  and  M.  Fink,  ‘Time  reversal  in  a  waveguide:  Study  of  the 
temporal  and  spatial  focusing,”/  Acoust.  Soc.  Am.,  vol.  107,  pp.  2418- 
2429,  2000. 

[16]  A.  Derode,  A.  Toupin,  J.  de  Rosny,  M.  Tanter,  S.  Yon  and  M.  Fink, 
‘Taking  advantage  of  Multiple  Scattering  to  Communicate  with  Time- 
Reversal  Antennas,”  Phys.  Rev.  Letters ,  vol.  90,  2003. 

[17]  T.  Strohmer,  M.  Emami,  J.  Hansen,  G.  Papanicolaou  and  A.  Paulraj, 
‘Application  of  time-reversal  with  MMSE  equalizer  to  UWB  communi¬ 
cations,”  to  appear  in  Proceedings  of  Globecom  Conference  2004,  2004. 


18 


[18]  P.  Kyritsi,  G.  Papanicolaou,  P.  Eggers  and  A.  Oprea,  ‘MISO  time 
reversal  and  delay-spread  compression  for  FWA  channels  at  5  GHz,” 
IEEE  Antennas  and  Wireless  Propagat.  Lett „  vol.  3,  no.  6,  pp.  96-99, 
2004. 

[19]  Y.  Chang,  H.R.  Fetterman,  I.L.  New  berg  and  S.K.  Panaretos,  ‘Mi¬ 
crowave  phase  conjugation  using  antenna  arrays ?  IEEE  Trans.  Microwave 
Theory  and  Techniques ,  vol.  46,  no.  11,  pp.  1910-1919,  1998. 

[20]  J.  Tuovinen,  G.S.  Shiroma,  W.E.  Forsyth  and  W.A.  Shiroma,  ‘Multi- 
path  communications  using  a  phase-conjugate  array”  IEEE  MTTS  Int. 
Microwave  Symp.  Digest ,  vol.  3,  pp.  1681-1684,  2003. 

[21]  F.  Dietrich,  R.  Hunger,  M.  Joham  and  W.  Utschick,  ‘Linear  precoding 
over  time-vaiying  channels  in  TDD  systems,”  IEEE  International  Conf. 
on  Acoustics,  Speech  and  Signal  Processing  ICASSP  ’03 ,  vol.  5,  pp.  117- 
120,  2003. 

[22]  R.  Esmailzadeh  et  ah,  ‘Prerake  diversity  combining  in  time-division 
duplex  CDMA  mobile  communications”  IEEE  Trans.  Veh.  Techn vol. 
48,  no.  3,  pp.  795-801,  1999. 


[23]  A.  Scaglione,  P.  Stoica,  S.  Barbarossa,  G.  Giannakis  and  H.  Sampath, 
‘Optimal  designs  for  space-time  linear  precoders  and  decoders,”  IEEE 
Trans.  Signal  Proc .,  vol.  50,  no.  5,  pp.  1051-1064,  2002. 

[24]  G.  Forney  and  M.  Eyuboglu,  ‘Combined  equalization  and  coding  using 
precoding”  IEEE  Comm.  Mag.,  vol.  29,  no.  12,  pp.  25-34,  1991. 

[25]  W.  Hirt  and  J.  Massey,  ‘Capacity  of  the  discrete-time  Gaussian  channel 
with  intersymbol  interference,”  IEEE  Trans.  Info.  Theory,  vol.  34,  no.  3, 
pp.  380-388,  1998. 

[26]  D.  Porcino  and  W.  Hirt,  ‘Ultra-wideband  radio  technology:  potential 
and  challenges  ahead,”  IEEE  Commun.  Mag.,  vol.  41,  no.  7,  pp.  66-74, 
2003. 

[27]  C.  Oestges,  A.D.  Kim,  G.  Papanicolaeou,  A.J.  Paulraj,  ‘Characterization 
of  space-time  focusing  in  time-reversed  random  fields,”  IEEE  Trans. 
Antennas  and  Propagat ,  vol.  53,  no.  2,  2005  (in  press). 

[28]  S.M.  Em  ami,  J.  Hansen,  A.D.  Kim,  G.  Papanicolaeou,  A.J.  Paulraj, 
D.  Cheung,  C.  Prettie,  ‘Predicted  time  reversal  performance  in  Wireless 
Communications  using  channel  measurements,”  IEEE  Comm .  Letters  (m 
press). 

[29]  M.  Emami,  M.  Vu,  J.  Hansen,  G.  Papanicolaou,  and  A.  Paulraj, 
‘Matched  fi  Itering  with  rate  back-off  for  low  complexity  communications 
in  very  large  delay  spread  channels,”  in  Proceedings  of  the  Asilomar 
Conference  2004,  2004,  submitted. 

[30]  W.C.Y.  Lee,  Mobile  Communications  Engineering,  McGraw-Hill,  1998. 

[31]  A.  Paulraj,  R.  Nabar  and  D.  Gore,  Introduction  to  Space-Time  Wireless 
Communications ,  Cambridge  University  Press,  2003. 

[32]  D.A.  McNamara,  C.W.I.  Pistorius  and  J.A.G.  Malherbe,  Introduction  to 
the  Uniform  Theory  of  Diffraction,  Artech  House,  1990. 

[33]  A.  Ishimaru,  Wave  propagation  and  scattering  in  random  media  (vol 
2),  Academic  Press,  1978. 

[34]  P.  Blomgren,  G.  Papanicolaou  and  H.  Zhao,  ‘Super-resolution  in  time- 
reversal  acoustics,”/  Acoust.  Soc.  Am.,  vol.  Ill,  pp.  230-248,  2002. 

[35]  J.  Hansen  and  A.  Paulraj,  Ttesign  approach  for  a  time  reversal  test  bed 
for  radio  channels,”  in  Proceedings  of  the  Eusipco  2004,  2004,  accepted. 

[36]  IEEE  802.16  Broadband  Wireless  Access  Working  Group,  ‘The  SUI 
channel  models,”  Tech.  Rep.  IEEE  802.1 6.3 c-01/29r5,  2003. 


IEEE  COMLET,  VOL.  XX  NO.  Y,  MONTH  2002 


l 


Predicted  Time  Reversal  Performance  in  Wireless 
Communications  Using  Channel  Measurements 

S.  M.  Emami,  Student  Member,  IEEE ,  J.  Hansen,  Student  Member,  IEEE ,  A.  D.  Kim,  G.  Papanicolaou, 
A.  J.  Paulraj,  Fellow,  IEEE ,  D.  Cheung,  and  C.  Prettie 


Abstract 

Using  broadband  radio  wireless  measurements  in  an  indoor 
environment  we  demonstrate  the  remarkable  space-time  focus¬ 
ing  properties  of  signal  transmission  with  time  reversal. 

I.  INTRODUCTION 

In  a  basic  time  reversal  experiment,  the  intended  receiver 
first  broadcasts  a  short  pilot  pulse.  The  transmitter  estimates 
the  channel  impulse  response  and  sends  the  time  reversed 
version  of  it  back  into  the  channel.  These  “time  reversed” 
waves  back  propagate  in  the  channel  by  retracing  their  paths 
and  focus  in  space  and  time  at  the  source,  the  intended 
receiver. 

In  a  channel  with  rich  scattering,  multipath  or  multiple 
scattering  is  exploited  by  TR  [1]  to  focus  broadband  signals 
tightly  in  space  and  time.  Moreover,  the  effective  channel 
obtained  through  TR  is  hardened  or  statistically  stable  [2], 
Spatial  focusing  means  that  the  spatial  profile  of  the  power 
peaks  at  the  intended  receiver  and  decays  rapidly  away  from 
the  receiver.  Temporal  focusing  means  that  the  channel  im¬ 
pulse  response  at  the  receiver  has  a  very  short  effective  length. 
Channel  hardening  is  the  phenomenon  of  tightening  of  the 
distribution  function  of  the  effective  channel  impulse  response. 
TR  works  well  in  systems  where  the  delay  spread  of  the 
channel  times  the  bandwidth  of  the  system  is  large  [3],  [2].  In 
this  paper,  we  demonstrate  TR  with  an  Ultra- Wideband  system 
in  an  indoor  environment. 

The  use  of  TR  has  three  main  benefits  in  communications. 
Temporal  focusing  significantly  shortens  the  effective  length  of 
the  channel.  For  example,  the  complexity  of  a  MLSE  equalizer 
is  exponential  in  the  length  of  the  channel.  TR,  thus,  reduces 
the  complexity  of  the  equalization  task. 

A  more  important  advantage  of  the  TR  technique  is  spatial 
focusing  [2].  Spatial  focusing  results  in  very  low  co-channel 
interference  in  a  multi-cell  system.  This  results  in  a  very 
efficient  use  of  bandwidth  in  the  overall  system. 

Another  advantage  of  using  TR  is  the  hardening  of  the 
effective  channel,  which  means  that  TR  results  in  a  high 
diversity  gain.  The  statistics  of  the  time  reversed  channel 
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are  different  from  the  actual  channel.  Specifically  the  time 
reversed  channel  has  a  much  smaller  variance  than  the  physical 
channel  itself  [2], 

Fink  and  his  coworkers  [1],  [3]  (and  references  therein) 
have  conducted  an  extensive  number  of  experiments  with  TR 
in  ultrasound  and  have  shown  large  gains  in  spatial  focusing 
in  the  presence  of  rich  scattering.  Also,  TR  methods  have 
been  used  for  underwater  acoustic  and  ultrasound  communi¬ 
cations  [4],  [5],  [6],  [3].  TR  experiments  using  radio  waves 
rather  than  acoustic  waves  have  not  yet  been  done.  Here, 
we  shall  perform  an  initial  demonstration  of  the  effects  of 
TR  methods  applied  to  radio  wireless  measurements  taken 
from  an  indoor  environment  in  the  band  2-8GHz.  Based  on 
these  measurements,  we  predict  the  performance  of  TR  in 
such  a  setting.  We  demonstrate  temporal  and  spatial  focusing 
for  a  SISO  link  and  evaluate  channel  hardening  in  terms  of 
the  empirical  distribution  of  the  received  signal  strength.  Our 
prediction  of  the  TR  performance  is  exact  if  there  is  no  channel 
measurement  error  at  the  transmitter  and  any  imperfections  in 
the  transmit  and  the  receive  chain  of  the  system  are  not  too 
severe. 

In  Section  II  we  explain  how  data  from  static  channel 
measurements  can  be  used  to  predict  the  performance  of  a 
TR  system,  and  we  define  suitable  performance  measures.  In 
Section  HI  we  describe  how  the  measurements  were  made.  In 
Section  IV  we  discuss  the  results  and  finally  conclude  with 
Section  V. 

II.  TR  Formulation 

Consider  a  transmitter-receiver  pair.  In  TR  the  transmitter 
uses  the  time  reversed  complex  conjugate  of  the  channel 
impulse  response  as  the  transmit  prefilter.  Denote  the  channel 
impulse  response  by  h( ro,  r),  where  r0  is  the  receiver  location 
and  r  is  the  delay  variable.  If  the  transmitter  uses  /i*(ro,  — t) 
as  the  prefilter,  the  effective  channel  to  any  location  r  is  thus 
given  by  the  time  reversed  field : 

s(r,  r)  =  h*(r0,  -r)  *  h(r,  r)  (1) 

where  *  denotes  convolution  with  respect  to  delay.  Note  that 
a  convolution  with  a  time-reversed  signal  is  equivalent  to  a 
correlation. 

In  order  to  study  the  effects  of  TR  we  measure  the  channel 
impulse  response  between  the  transmitter  and  a  receiver.  We 
repeat  the  measurement  by  holding  the  transmitter  fixed  and 
changing  the  position  of  the  receiver  over  a  square  grid.  We 
take  one  comer  of  the  grid  as  the  reference  receiver  (the 
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intended  receiver)  and  use  its  impulse  response  as  the  transmit 
prefilter  h*( r0,  —  r).  The  time  reversed  field  at  any  point  on 
the  square  grid  can  now  be  computed  based  on  the  channel 
measurements  and  by  using  Eqn.  1.  This  is  equivalent  to 
having  the  transmitter  perform  the  prefiltering  with  the  TR 
filter. 

The  power  of  the  signal  as  a  function  of  both  space  and 
delay  is  computed  and  investigated.  As  a  metric  for  spatial 
focusing  we  define  the  following  quantity: 

«(r)  =  max  |s(r,  r)|2  (2) 

This  quantity  represents  the  power  of  the  strongest  tap  at  a 
receiver  located  at  r.  We  expect  «(r)  to  peak  at  tq  and  to 
decay  rapidly  with  increasing  distance  from  the  origin. 


III.  Measurement  Setup 

Measurements  were  conducted  by  Intel  Corporation  at  off- 
peak  hours  to  ensure  channel  stationarity.  The  environment  is 
an  office  space  (40m  x  60m)  with  many  cubicles.  Measure¬ 
ments  span  the  bandwidth  2-8  GHz  with  3.75  MHz  frequency 
resolution.  From  the  data  we  estimate  that  the  coherence 
bandwidth  of  the  channel  is  20  MHz  Antennas  are  vertically 
polarized.  The  virtual  grid  on  which  the  receiver  is  moved  has 
a  distance  of  Aq/4  where  Aq  is  the  wave  length  of  the  mid 
frequency  of  the  measurements  (5  GHz). 

The  receive  antenna  is  moved  to  a  different  location  with 
a  precise  robotic  positioner.  At  each  antenna  position,  the 
channel  is  measured  using  a  vector  network  analyser.  The 
measurements  are  corrected  to  compensate  for  the  system 
components  (including  cable,  gain  stages,  and  antennas).  The 
height  of  the  transmit  antenna  is  about  2.5m  and  that  of  the 
receive  antenna  is  lm  above  the  floor.  The  channel  impulse 
response  was  computed  by  taking  the  inverse  DFT  of  mea¬ 
surement  data. 

IV.  Experimental  Results 

We  evaluated  «(r)  for  11  different  scenarios.  Here  we 
present  two  cases.  Case  I  is  an  essentially  line-of-sight  scenario 
with  a  separation  of  3m.  Case  II  has  a  separation  of  1  lm  and 
is  a  typical  non-line-of-sight  situation.  The  results  for  the  two 
cases  discussed  here  have  the  worst  focusing  and  yet  they  are 
very  similar  in  their  outcomes  with  all  the  other  scenarios. 

In  the  3-D  figures,  the  square  grid  spans  a  region  of 
7Aq  x  7A0.  A  3D  plot  of  «(r)  of  the  two  cases  is  shown 
in  Fig.  1.  We  see  that  spatial  focusing  works  fine  in  both 
scenarios.  In  neither  of  them,  however,  the  peak  is  isotropic; 
both  peaks  have  one  direction  in  which  they  Ml  off  faster, 
and  another,  in  which  the  decay  is  slower.  The  structure  of 
the  peak  carries  some  information  about  the  geometry  of  the 
environment,  i.e.,  about  directions  which  show  faster  and  those 
which  show  slower  decorrelation  in  space. 

We  also  observe  that  the  signal  power  level  is  at  least  10 dB 
lower  at  a  distance  of  7A0  than  its  value  at  the  receiver. 
This  demonstrates  the  spatial  decorrelation  very  well.  We  can 
conclude  that  channel  impulse  responses  show  inherent  quasi¬ 
orthogonality  in  space,  which  can  be  used  for  interference 
suppression.  In  Fig.  2  we  demonstrate  the  compression  of  the 


(a)  Case  I 


Distance  in  wavelength  Dtetance  10  '“"3th 

(b)  Case  II 

Fig.  1.  Spatial  Focusing  «(r).  One  shot  spatial  field  realizations  for  a)  the 
line-of-sight  b)  the  non-line-of-sight  scenario. 


received  pulse  in  the  time  domain.  Displayed  are  the  channel 
impulse  response  magnitude  for  the  LOS  scenario  (a),  the 
respective  time-reversed,  i.e.,  compressed  impulse  response 
magnitude  (b),  the  impulse  response  for  the  NLOS  scenario 
(c),  and  the  compressed  one  for  this  case  (d).  The  effective 
length  of  the  channel  is  significantly  shortened.  Still,  for 
the  LOS-case,  comparably  strong  side  lobes  remain  visible. 
The  impulse  response  of  the  NLOS  channel  is  much  better 
compressed  in  time.  There  is  some  intersymbol  interference  in 
the  time  compressed  field  and  one  would  need  to  equalize  the 
received  signal.  However,  the  time  compressed  channel  has  a 
significantly  shorter  effective  length  and  is  thus  less  expensive 
to  equalize. 

In  Fig.  3  we  investigate  the  hardening  of  the  channel.  We 
consider  the  distribution  of  the  strongest  tap  of  the  channel 
impulse  response  and  of  time-reversed  field.  The  dashed  curve 
represents  the  distribution  of  the  zero-th  tap  of  the  channel 
impulse  response,  and  the  solid  curve  represents  that  of  the 
compressed  signal.  The  figure  illustrates  that  the  time  reversed 
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(a)  Case  I 


(b)  Case  I 


(c)  Case  II 


(d)  Case  II 


Fig.  3.  Distribution  of  s(r.O)  (solid)  and  that  of  /i(r,0) 

field  is  much  less  random  (i.e,,  has  more  diversity)  than  the 
channel  response  itself. 

V.  Conclusions 

Using  the  indoor  measurements  discussed  above,  TR  results 
in  strong  spatial  focusing  and  time  compression.  At  about  7 
wavelengths  from  the  target  the  signal  power  reduces  by  at 
least  lOdB  from  its  value  at  the  target.  The  channel  impulse 
response  is  well  shortened  at  the  target  as  well.  Channel 
hardening  or  diversity  gain  is  also  observed.  In  a  practical 
setting  the  results  may  degrade  if  the  transmitter  does  not 
perfectly  know  the  channel  or  if  the  transmitter  circuitry  has 
nonlinearities.  Our  calculations  with  the  indoor  radio  wireless 
measurements  obtained  by  Intel  suggest  that  the  benefits  of 
TR,  which  have  been  observed  in  ultrasound  and  underwater 
sound  TR  experiments,  could  also  be  gained  in  wireless 
communications. 
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Fig.  2.  Magnitude  of  channel  impulse  responses  (a  and  c)  and  the  time 
compressed  impulse  response  (b  and  d) 


